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ABSTRACT 

Aims. We are studying the dust properties of four low metallicity galaxies by modelling their spectral energy distributions. This 
modelling enables us to constrain the dust properties such as the mass, the temperature or the composition to characterise the global 
ISM properties in dwarf galaxies. 

Methods. We present 870 /um images of four low metallicity galaxies (NGC 1705, Haro 11, Mrk 1089 and UM 311) observed with 
the Large APEX BOlometer CAmera (LABOCA) on the Atacama Pathfinder Experiment (APEX) telescope. We model their spectral 
energy distributions combining the submm observations of LABOCA, 2MASS, IRAS, Spitzer photometric data and the IRS data for 
Haro 11. 

Results. We find that the PAH mass abundance is very low in these galaxies, 5 to 50 times lower than the PAH mass fraction of 
our Galaxy. We also find that a significant mass of dust is revealed when using submm constraints compared to that measured with 
only mid-IR to far-IR observations extending only to 160 /im. For NGC 1705 and Haro 11, an excess in submillimeter wavelengths 
is detected when we use our standard dust SED model. We rerun our SED procedure adding a cold dust component (10 K) to better 
describe the high 870 jim flux derived from LABOCA observations, which significantly improves the fit. We find that at least 70 
% of the dust mass of these two galaxies can reside in a cold dust component. We also show that the subsequent dust-to-gas mass 
ratios, considering HI and CO observations, can be strikingly high for Haro 1 1 in comparison with what is usually expected for these 
low-metallicity environments. Furthermore, we derive the star formation rate of our galaxies and compare them to the Schmidt law. 
Haro 1 1 falls anomalously far from the Schmidt relation. These results may suggest that a reservoir of hidden gas could be present 
in molecular form not traced by the current CO observations. While there can be a significant cold dust mass found in Haro 11, the 
SED peaks at exceptionally short wavelengths (36 /um), also highlighting the importance of the much warmer dust component heated 
by the massive star clusters in Haro 11. We also derive the total IR luminosities derived from our models and compare them with 
relations that derive this luminosity from Spitzer bands. We find that the Draine & Li (2007) formula compares well to our direct IR 
determinations. 

Key words. galaxies:ISM - galaxies: dwarf - Infrared:ISM - ISM:dust,extinction 



1. Introduction 

The understanding of the evolution of a galaxy requires knowl- 
edge of the roles of the different actors controlling the evolution 
of the Interstellar Medium (ISM) and the subsequent feedback 
on star formation activity. Despite its low fraction of the total 
mass of a galaxy (less than 1%), dust plays a prominent role in 
the heating and cooling of the ISM and thus tightly influences 
the overall physics of a galaxy. Since dust absorbs the stellar ra- 
diation and reemits it in a wide range of wavelengths, the star 
formation rate (SFR) as well as other fundamental parameters 
of a galaxy, such as its age, can be indirectly studied through the 
dust emission itself. The Spectral Energy Distribution (SED) of a 
galaxy is its spectral footprint from which we can study the phys- 
ical processes taking place in the galaxy since it synthesises the 



contribution of all its components to the emission of the galaxy. 
Using this tool, we can peer into the window of the integrated 
history of the galaxy and disentangle the various physical actors 
(stars, HII regions, molecular clouds) and processes (stellar ra- 
diation, dus t emission) involved (Drain e~et al.| |2007; Gall iano | 
|et al.||2008l see also § 5 of this paper). 

While dust hinders the interpretation of ultraviolet (UV) and 
optical wavelengths, in the Mid Infrared (MIR), Far Infrared 
(FIR) and submillimetre (submm) wavelengths, dust emission 
and absorption properties expose different physical environ- 
ments, from the most vigorous star form ation and AGN activ- 
ity (e.g. |Gordon||1995[|Wu et al.||2007|) to the more q uiescent 



diffuse media ( [Bernard et al.||1996[|Arendt et al.| [T998]). Many 



processes linked to star formation such as stellar winds ( Hoefner 



2009 ), supernovae shocks, photodestruction by high-mass stars 
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etc. can also affect the spatial distribution and the local properties 
and abundance of the different dust components of a galaxy such 
as Poly cyclic Aromatic Hydrocarbons (PAHs, O'Halloran et al.| 



flLisenfeld et al |[200l]|Bottner et alj[200^|Dumke eTaLl|2004 



2006), amorphous carbon grains, silicates or composite grains, 



manifesting themselves in the MIR to submm wavelengths. 

Studying the interplay between galaxy properties and metal 
enrichment is crucial to understand galaxy evolution. The metal- 
licity of a galaxy is deeply linked with the dust properties of 
the ISM and its substructures such as HII regions and molec- 
ular clouds, but just how it affects the ISM is currently poorly 
known. Dwarf galaxies in the Local Universe, are metal-poor 
galaxies, and are thus convenient laboratories to study the effects 
of metallicity (Z) on the gas and dust. They exhibit a wide va- 
riety of physical conditions, and their star formation properties 
and ISM represent the closest analogs to proto-galaxies of the 
early universe. Indeed, dwarf galaxies are small and may com- 
pared t o high redshift galaxies w hich also present lower metal- 
licities ( |Lara-L6pez et al.| |2009). They are also considered to be 
the building blocks of much larger and more metal-rich galaxies 
(Review by Tosi, 2003 ). They also show analogies with Gamma 
Ray Bursts (GRB) hosts whose ISM usually exhibit moderate 



chemical enrichment with a median metallicity of 1/10 Z G ( |Chen 
[eTaL|[2009] ). 

They finally show evidence for older stellar populations than 
their metallicity suggests (e.g. |Aloisi et al.|[T99 8), posing enig- 
matic issues for galaxy evolution models. Many studies have 



been carried out to grasp this apparent paradox. Lisenfeld & 
|Ferra ra (1998) confirmed that the dependence of the dust-to-gas 
mass ratio (D/G) in low metallicity galaxies was a function of 
metallicity using IRAS observat ions.! James et al.| ( [2002] ), [Walter 



|et al.| ( [2007] ) and |Hirashita et al. ( |2008| ) concluded likewise using 



JCM77SCUBA submm, Spitzer MIR/FIR and AKARI (FIR) ob- 
servations. Finally, Galliano et al. (2008) observed some system- 
atic deviations between dust abundances of very low metallic- 
ity systems and what is expected for supernova-condensed dust. 
At MIR wavelengths, low metallicity systems also show promi- 
nent differences in the dust properties compared to the more 
metal-rich systems. For example, PAH features are strikingly di- 
minished as metallicities drop (e.g. |Madden| |2005[ [Engelbracht 
et~aLl [20051 |Wu et al.[ [20071 |Engelbracht et al.f |2008) com- 



pared to metal-rich galaxies, in spite of th e role the smalles t 
grains play in the energy balance of galaxies ( Rubin et al. 2009| ). 
Some studies suggest that PAH emission depends on the hard- 



ness or strengh of the illuminating radiation field ( Madde n et al. 
2006] [EngeTbracht et al.||2008l|Gordonet^l|2008HBendo et al. 
2008J). The consequence of lowering the metallicity of a galaxy 
is the decrease in dust opacity resulting in harder and stronger 
radiation fields. The dearth of PAHs in low metallicity galaxies 
has also been explained by the d estructive effects of supernovae 
2008|) or by the delayed injection of 



(O'Halloran et al. 2006 



PAHs by AGB stars ([GaT 



ianoetal. 2008). 



Broad wavelength coverage of the MIR to submm regime 
is imperative to constrain the modelling of the observed SEDs, 
leading to a better comprehension of the dust properties of galax- 
ies. Since Spitzer only observes dust emission at wavelengths 
shorter than 160 pm, submm data are necessary not only to 
enlarge the wavelength coverage at longer wavelengths to ver- 
ify the dust models but also because the potential reservoir of 
cold grains (< 15K), which contribute to this submm flux, may 
account for a significant amount of mass. Only a handful of 
galaxies of the Local Universe have been studied using submm 
ground-based instruments (e.g. /CM77SCUBA). When submm 
observations of dwarf or late-type galaxies are studied, an ex- 
cess in the dust SEDs is often found in the mm/submm domain 



Gallia noeFal.1 |2003[ |2005[ |Marleau end! |2006[ |Bendo et al. 



2006). This excess can be interpreted as very cold dust (<10K), 
in which case more than 50% of the total dust mass of these 
galaxies should reside in a very cold component. Cold dust is 
also needed to explain the break in the Gas to Dust mas ratio 
as a function of metallicity relaion for low-metall icity galax- 
ies ( [Galliano et"aLl|2008[[Munoz-Mateos et al.[|2009| ). The pres- 
ence of this cold dust component is still a contentious issue in 
the ISM community and will have important consequences on 
our comprehension of ISM properties of low metallicity environ- 
ments. |Lisenfeld etaL](|2001| ), IReach et al |(|1995| ), |Dumke et aL] 
(2004|), |Bendo et al.| ( |2006| ) or |Meny et al.| ( |2007] ) suggested that 
changes in dust emission properties (changes in dust emissivity 
or resonances related to dust impurities) should be responsible 
for boosting submm emission above the 15-20K thermal emis- 
sion expected at these wavelengths. However, not all low metal- 
licity galaxies show submm excess, as shown recently by the 
observations of the nearby Local Group Galaxy IC10 (Parkin et 
al. 2009 in prep), where the main two star forming regions were 
isolated with ISO, Spitzer and 850 pm observations, the SEDs 
were modeled without invoking a very cold dust component. 
Moreover, Draine et al. (2007 ) showed that their observations 
of mostly metal-rich galaxies can largely be reproduced by dust 
models which do not account for a very cold dust component, 
even in their limited number of cases where submm observa- 
tions are present. Studies using submm observations for a wider 
range of metallicity values are necessary to check the relevance 
of these conclusions for low metallicity environments. 

In this paper, we present the first APEX/ LABOCA 870 pm 
observations of dwarf galaxies: 1 extended galaxy (NGC 1705) 
and 3 compact sources (Haro 11, UM 311, Mrk 1089) (see 
Table [2]). We have combined these data with Spitzer and/or IRAS 
observations to produce global SEDs that we use to model the 
dust properties. The sample is small but covers a wide range 
of metallicities, from ~l/9 Z for Haro 11 to ~l/3 Z© for 
NGC 1705. It also presents varied morphologies, size scales and 
characteristics: resolved or compact galaxies, disturbed and even 
merging environments. 

We describe the sample in § 2 and the observations and data 
reduction in § 3 and the images and photometry in § 4. In § 5, 
we present the SED modelling and discuss the results in § 6. 



2. The sample 

The four galaxies studied in this paper were chosen because of 
their diversity in morphology, distance, metallicity and star for- 
mation activity: 

NGC 1705 - This is a well studied Local Group galaxy 
whi ch is in the Spitzer Infrared Nearby Galaxie s Survey (SINGS; 
|Kennicutt et"aT] |2003] ). |Heckman & Leitherer| ([1997]) found the 
luminosity of the galaxy in the UV to be dominated by a cen- 
tral bright 10 5 M o Super Star Cluster (SSC) also at the origin 
of a galactic outflow. This SS C shows similar p roperties to the 
distant gamma ray burst hosts ( [Chen et aL) 2007| ). The Ha emis- 
sion extends over the entire optical emission of the galaxy ([Gil 



de P az et al.[ |2003| ) while the HI e mission is exc eptionally ex- 
tended beyond the optical emission (Meurer et al.||1998"]) and lie 



on either side of the SSC. Two off-nuclear regions called Dl and 
D2 ( Fig.[T]) can be s een to dominate the MIR and FIR dust emis- 
sion ( Cannon et al.||2006| ). Spitzer IRS spectroscopy reveals the 
PAH emission ori ginating toward regio n Dl but not toward the 
SSC or region D2 ( [Cannon et aL) |20U6] ). 
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Name 


Ra (2000) 


Dec (2000) 


12+log(0/H) 


Distance (Mpc) 


M(HI) 


Apparent size 


NGC 1705 
Haro 11 

Mrk 1089 
UM311 


04h54' 13.5" 
OOh 36' 52.5" 
05h01' 37.8" 
Olh 15' 34" 


-53° 21' 40" 
-33° 33' 19" 
-04° 15' 28" 
-00° 51' 32" 


8.46 (a) 
7.9 {b) 
8X) & 
8.3 ^ 


4.7 (a) 
92 (c) 
59.8 {d) 
21.3 (J) 


5.1 x 10 7 (1) 

~ 10 8 (2) 
2.7 xlO 10 ® 
2.3 x 10 9 (4) 


1.9' x 1.4' 
0.5' x0.5' 
0.61' xO.23' 
0.11' xO.ll' 





Meurer et al. 



(1992), b) Bergvall & Qstlin (2002 ), c) Bergvall et al. (2000), d) Hopkins 



Table 1. General properties of the sample: a) 

|etaLl ( [2QQ2l » 

The HI masses were scaled according to the distances we used. Re ferences for HI dust mass : 1) 8.9 x 10 7 M© for D=6.2 Mpc 



1998 ), 2) Upper limit derived from the H I detection (B ergvall & 6stlin[ |2002| ), 3) 2.1 x 10 10 M Q for D=53 Mpc 



(Meu rer et al . A1 

(Williams et al.||1991|), 4) 3.1 x 10 9 M Q for D=24.7 Mpc ([Smoker et aT.||2000|). 



Haro 11 - Also known as ESQ 350-IG038, this is the most 



distant galaxy of the sample (92Mpc; Bergv all efaLl |2000| ) 
It posses ses characteristics of a n extreme starburst with h FIR 



10 11 L© ( Sanders et al. 20031), making it a luminous infrared 
galaxy (LIRG) with a high star formation rate of ~ 25 M Q yr" 1 as 
determined from Ha, radio continuum, FIR and hard X-ray ob- 



servations ( Grimes et al. 2007|). Broadband images of Ha show 



three bright star-forming condensations with unrelaxed kine- 
matic structure and faint extended shell structures in the outer re- 
gions of the galaxy, all suggesting an ongoing merger ( Bergvall 



& Ostlin 



2002; Ostlin et al. 



[1999]). Haro 11 is a moderately 
strong radio source (essentially free-free continuum) with ex- 
tended continuum emission at 6 and 20 cm (Heisler & Vaderl 



1995 ). It is a very metal poor galaxy (Z-l/7 Z Q ) that seems to 
have a very little neutral hydrogen, an unusually high ratio be- 
tween blue luminosity and HI mass and little observed molecular 
gas ( [Bergvall et aLl[2000| ). 

Mrk 1089 - Mrk 1089 is a Wolf-Rayet (WR) galaxy ( |Kunth| 
& Schild 1986| ) and is the most luminous of the eight mem- 
bers of the Hickson Group 31 (Hickson, 1982). The morphol- 
ogy of the group is very disturbed with tidal interactions be- 
tween Mrk 1089 (HCG31 C) and the galaxy NGC 1741 (HCG31 
A)(Fig. [2J:). The two ga laxies present simi lar kinematics sug- 
gesting a single entity (Ri cher et al.| |2003| ). At the location of 
their interaction is a very strong 24 fim source and high levels 
of infrared emission c an be found throughout the whole group 
( [Johnson et aL| [2007] ). Mrkl089 and NGC1741 are referred to 
collectively as NGC 1741 in the catalog of WR galaxies of |Conti 
(|1991|). Nevertherless, many surveys ([Rubin et al. 1990 



Conti 



etal.l |1996; Igl esias-Paramo & Vilchez||l 997 [| Johnson & Conti 



2000 ) showed that the interacting system was undergoing a star- 
burst which was attributed to HCG31 C. We thus decide, in 
this paper, to designate the interacting system compounded of 

HCG31 A and HCG31 C together as Mr k 1089. 

T99T> 



UM 311 - This compact HII galaxy (Terlevich et al. 



is located between the pair of spiral galaxies NGC 450 and 
UGC 807, 2 galaxies which were once thought to be interact- 
ing but have now been demonstrated to be physically separated 
(Rubin & Ford[ |1983| ). There are three very bright sources of 
compact HII emission between the two galaxies, UM 311 being 
the brightest (Fig. [2]l). The galaxy has been misidentified as a 
projected galactic star due to its quasi- stellar and compact mor- 
phology. Its Ha luminosi ty and equivalent w idth are remarkably 
high for an HII galaxy ( Gus eva et al. 1998 ). 



3. Observations and data reduction 

3.1. LABOCA 

LABOCA is a multi-channel bolometer array for continuum 
observations at 870 /mi, built by MPIfR (Max-Planck-Institut 



fur Radioastronomie, Bonn, Germany) and mounted on APEX 
(Atacama Pathfinder Experiment), a 12-meter radio telescope 
of ESO, Onsala and MPIfR. The array consists of 295 channels. 
LABOCA on APEX has a total field of view of 11^4x1 YA and a 
full width half maximum (FWHM) of its point spread function 
(PSF)~18.2". 

About 30 hours of observations were taken from November 
9th 2007 to November 20th 2007. The four galaxies are smaller 
than the FOV of LABOCA. Basic spiral patterns with 4 point- 
ings were combined with a raster mapping mode (raster- spiral) 
to completely fill the array (homogenous sampling of each map), 
to obtain enough off-source position for background substrac- 
tion and to effectively execute long integration times (8.7h, 8.3h, 
4h and 7.6h for Haro 11, UM 311, Mrk 1089 and NGC 1705 
respectively). Flux calibration was performed through the obser- 
vations of the planets Uranus and Mars and the sources HLTAU, 
J0423-013, J0050-095, J0006-064, V883-ORI, PMNJ0403- 
3605, PMNJ0106-4034, PKS0537-441, N2071IR. We used the 
BOA package (BOlometer Array Analysis Software; Schuller 
et al., in prep) to reduce the data. The sotware was devel- 
oped through a collaboration of scientists from the MPIfR, 
AlfA (Argelander-Institut fur Astronomie, Bonn, Germany), 
AIRUB (Astronomisches Institut der Ruhr-Universitat, Bochum, 
Germany), and IAS (Institut d'Astrophysique Spatiale, Orsay, 
France). 

The atmospheric attenuation was determined via skydips 
every hour. Opacities at 638 /im range from 0.103 to 0.353 
throughout the observing run. Flat-fielding normalisation is ap- 
plied by dividing the signals by the bolometre gains supplied 
by the observatory and thus removing the bolometre gain varia- 
tions. From the 18th of October up to the 15th of November, 33 
LABOCA pixels were shadowed by a plate in the beam. These 
pixels are removed. We also masked out stationary points and 
data taken outside reasonable telescope scanning velocity and 
acceleration limits as well as dead or noisy channels. During the 
commissioning period of LABOCA, correlated noise was found 
between groups of channels sharing some parts of the electronics 
(amplifier boxes or cables). As our data are not all point sources, 
it is necessary to carefully differentiate correlated noise from ex- 
tended emission. This correlated noise is subtracted from each 
map. Moreover, we performed a correction to suppress spikes 
and flatten the 1/f noise of the FFT to remove the noise caused 
by thermal variations. Finally, each reduced scan was gridded 
into a weighting map. These weights were built by calculating 
the rms of each time series which contributes to a given region 
of the map. 

The major sources of uncertainty are the calibration uncer- 
tainty and uncertainty resulting from the background variation. 
We estimate the average uncertainty levels to be of ~ 15 %. All 
of the galaxies were detected with a root mean square (rms) esti- 
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NGC1705 (Hot) b) NGC1705 (8 pm + Ha contours) 




4 h 54 m 19 s 4 h 54 m 17 s 4 h 54 m 14 s 4 h 54 m 12 s 4 h 54 m 10 s 4 h 54 m 7 s 
Right Ascension (J2000) 



D1 SSC D2 




4 h 54 m 19 s 4 h 54 m 17 s 4 h 54 m 14 s 4 h 54 m 12 s 4 h 54 m 10 s 4 h 54 m 7 s 
Right Ascension (J2000) 



c) 



NGC1705 (HI) 



d) 



§ -53°2 




NGC1705 (8 pm + HI contours) 




4 h 54 m 2 gs 4 h 54 m 24 s 4 h 54 m 19 s ^m^s ^g^gs 

Right Ascension (J2000) 



4 h 54 m 29 E 4 h 54 m 24 s 4 h 54 m 1 9 s 4 h 54 m 1 4 s 4 h 54 m 1 s 4 h 54 m 5 s 4 h 54 m 
Right Ascension (J2000) 



Fig. 1. a) Ha image of NGC 1705. Note that the SSC is masked in the image ( [Gil de Paz et al.||2003| ). b) NGC 1705 obs erved at 
8 pm with Spitzer. Ha contours are ove rlaid. Note the SSC and the two bright off-nuclear HII regions labeled Dl and D2 ( |Cannon| 
etaL|[2006| ). c) HI image of NGC 1705 ( |Meurer et aLl[2006| ). d) NGC 1705 observed at 8 /mi with Spitzer. HI contours are overlaid. 



mated to be lower than 10 mJy/beam for three of our galaxies but 
as high as 15 mJy/beam for Mrk 1089. Some anomalous pixels 
at the edges of the final maps are the results of poor coverage at 
the edges of the map. The anomalous pixels are ignored in the 
analysis. The 870 pm LABOCA images are presented in Fig. [2] 



Request (AOR) keys of each observation. IRAC fluxes are un- 
certain at the ~ 10% level due to systematic effects. The data 
we obtained were c ompared and verified with those publish ed in 
IJohnson et al.1j2007) for Mrk 1089 and |Cannon et al.] ( |2006] ) for 
NCG 1705. Images are presented in Fig. [2] 



3.2. Spitzer data 



3.2.2. MIPS 



3.2.1. IRAC 

The IRAC bands cover 3.6, 4.5, 5.8, and 8 pm with a FWHM 
of the PSFs of 1.8, 1.8, 2.0 and 2.2 arcsec respectively 
(http://ssc.spitzer.caltech.edu/documents/SOM/). Two bands are 
imaged in pairs (3.6 and 5.8 microns; 4.5 and 8.0 microns). 
IRAC provides 5.2' x 5.2' images (pixel size of 1.2" x 1.2"). Our 
sources were observed with IRAC in dithering imaging mode 
and Mrk 1089, UM 311 and Haro 11 were obtained through the 
Spitzer data archive (post-basic calibrated data). For NGC 1705, 
the IRAC data were obtained throught the SINGS data deliv- 
ery page 1 . Table [2] summarizes the Astronomical Observation 



1 Details on the data reduction of the SINGS Fifth Enhanced Data 
Release can be found at http://data.spitzer.caltech.edu/popular/sings/ 



The MIPS bands cover 24, 70 and 160 pm with a 
FWHM of the PSFs of 6, 18 and 4 0" respectively 
(http://ssc.spitzer.caltech.edu/documents/SOM/). The galaxies 
were observed in MIPS scan mapping mode for NGC 1705 and 
MIPS Photometry/Super-Res mode for Haro 11, Mrk 1089 and 
UM311. 

All of the MIPS images (see Table [2] for AOR keys) were 
obtained through the Spitzer data archive (raw data) and reduced 
using the MIPS Data Analysis Tools ( [Gordon etaL] |2QQ5| ), ver- 
sion 3.10 along with additional processing steps. 

The individual 24 pm frames were first processed through 
a droop correction (to remove an excess signal in each pixel) 
and were corrected for non-linearity in the ramps. The dark cur- 
rent was then subtracted. Scan-mirror-position dependent and 
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NGC 1705 (3.6 /jm) NGC 1705 (8 //m) 

' ' I 00" 



04 h 54 m 20 s 



15 s 1(T 
a J2000 




NGC 1705 (24 /im) 



NGC 1705 (870 /im) 



b) 




04 h 54 m 20 s 



a J2000 

Haro 11 (3.6 iim) 





"J20QG 

Haro 11 (8/im) 



o 33 



o 



00 37 00 36 55' 




Haro 11 (24 /im) 



o 33 



00 h 37 m 00 s 36 m 55 s 




Haro 11 (870 A/m) 



00 h 37 m 00 s 36 m 55 s 



Fig. 2. 3.6 fim (4.5 for UM 311), 8 jim , 24 jim and 870 yum images of our four galaxies: a) NGC 1705, b) Haro 11, c) Mrk 1089 
(In the paper, the interacting system of NGC 1741 and Mrk 1089 is designated as Mrk 1089), d) UM 311 (The galaxy is identified 
as source number 3). North is upwards and East is to the left. The FWHM of the PSF of LABOCA is shown by a white circle. 8 fim 
(4.5 /mi for Haroll) contours are overlaid on the LABOCA images: a) 0.5, 1, 2.5 MJy/sr - (b) 0.34, 0.77, 8.58 and 9 MJy/sr - c) 
2, 2.5, 4.5, 8 MJy/sr - d) 4.6, 5.61, 7.8 MJy/sr. For UM 311, the photometric aperture chosen (54") encompasses the flux of the 3 



sources at the edges of the two spiral galaxies and is marked by the dashed circle. These sources are numbered according to Moles 



et al.| ( |1994| ). For NGC 1705, three individual regions can be isolated: the SSC in the center and two off-nuclear HII regions on both 



sides of this SSC. 
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c) 



M. Galametz et al: The SED of dwarf galaxies using LABOCA 
Mrk 1089 (3.6 /urn) Mrk 1089 (8 //m) 



-04° 16' 30" 





Mrk 1089 
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NGC 1741 J 




• 






W ■ 

• 



-04 1J 16 r 30' 



05 h 01 m 40 s 




05 h 01 m 40 s 



Mrk 1089 (24 jim) 



Mrk 1089 (870 //m) 



o 15' 30' 



-04° 16' 30" - 




05 11 01 111 40- 




d) 



UM311 (4.5 /mi) 



UM311 (8 /mi) 



UGC 807 



NGC 450 



: ♦ 



-00° 53' 



35 s 30 s 
a J2000 




35 s 30 b 
a J2000 




Fig. 2. continued 
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Name AOR key Spitzer/IRAC AOR key Spitzer/MlPS 



NGC 1705 5536000 - 5535744 5549056 - 5549312 - 14812416 
Haroll 4326400 4344832 

Mrkl089 11250432 11247360 

UM311 10392576 12562176 



Table 2. AOR keys of the Spitzer/IRAC and Spitzer/MSPS observations. 



independent flats were created in each AOR and applied to the 
data. Detector pixels that had measured signals superior to 2500 
DNs -1 (data numbers per second) in any frame were masked 
out in the following three frames to avoid latent images in the 
data. Planes were fit to the zodiacal light emission in the back- 
ground regions, and third order polynomials were fit to the back- 
ground in each leg of each scan map. We also measured scan- 
mirror-position dependent residual offsets in the backgrounds. 
They were subtracted from the data. For NGC 1705, we cre- 
ated preliminary mosaics of the data to identify transient objects 
(e.g. asteroids). We also performed a robust statistical analysis in 
which the values of cospatial pixels from different frames were 
compared to each other and statistical outliers to remove cosmic 
rays. We masked out these objects. A final mosaic was made 
with data calibrated into astronomical units. The calibration fac- 
tor for the 24 /mi data was given in Engelbrach t et al.|(2007| ) as 
4.54 + 0.18 x 10 -2 MJy sr" 1 [MIPS instrumental unitpV" 

In the 70 and 160 /mi data processing, we first fit ramps to 
the readouts to derive slopes. Readout jumps, cosmic ray hits and 
dark current were removed. The stim flash frames taken by the 
instrument were used as responsivity corrections. An electronic 
nonlinearity correction and an illumination correction were ap- 
plied. Short term variations in the signals ('drifts') were removed 
from all 70 /mi and 160 /mi data. To combine the photometry 
and scan map data for NGC 1705 without dealing with prob- 
lems of background offsets, we subtracted the background from 
the photometry map at this stage. We corrected all pixels affected 
by cosmic rays (see 24 /mi data treatment for details). Then, fi- 
nal mosaics were built, residual backgrounds were subtracted 
and the data calibrat ed. The 70 /mi calibration factors given in 
Gordon et alT] ( |2007| > are 702 ± 35 MJy sr" 1 [MIPS instrumen- 
tal unit] -1 for coarse-scale imaging and 2894 ± 294 MJy sr" 1 
[MIPS instrumental unit] -1 for fine-scale imaging. The 160 /mi 
calibration factor is given by Stansberry et al. (2007 ) as 41.7 ± 5 
MJy sr -1 [MIPS instrumental unit] -1 . An additional 70 /mi non- 
linearity correction given as f70(true) = 0.581 X f70(measured) 113 
( [Dale et aL| |2007j > was applied where the surface brightness ex- 
ceeded 66 MJy sr -1 . MIPS flux uncertainties were estimated to 
be 10% in the MIPS 24 /mi band and 20% in the MIPS 70 and 
160 /mi bands. The fluxes ca lculated compared t o those pub- 
lished in [Johnson et al.| ( |2007| ) for Mrk 1089 and |Cannon et aT] 
([2006]) for NCG 1705. 



tion channel (SL, 5.2-14.5 /mi; /tyA/l~60-120) an d the two high 
resolution channels: short-high (SH, 9.9-19.6 /mi; A/AA-600) 
and long-high (LH, 18.7-37.2 /mi; A/AA-600). IRSClean was 
applied to the individual files to mask out hot pixels, and 
the frames were coaddd together. We then used the Spitzer 
IRS Custom Extraction software (Spice - 2.0.1 version) to ex- 
tract the spectrum. We subtracted background contributions es- 
timated from zodiacal models using the software SPOT (18.0.1 
version - http://ssc.spitzer.caltech.edu/propkit/spot/). We finally 
subtracted a background estimated using zodiacal models in 
SPOT. The final error on the photometry derived from peak-up 
acquisition windows is about 15%, the value recommended by 
the IRS data handbook. Finally , the spec trum is shifted from a 
redshift of z=0.021 (Bergval l et al.| [2000) to rest wavelengths. 
The position of the slits of the IRS high resolution observations 
is shown in the panel on the left in Fig. [3] The spectrum with la- 
bels identifying the main emission features is shown in the panel 
on the right in Fig. [3] The 6.3, 7.6, 8.2 and 11.3 /mi PAH fea- 
tures are clearly detected as well as significant emission lines. 
Note that Nelll/Nell ratio=3.2 (calculated on the IRS spectrum), 
which is an indication of the hard radiation field dominated by a 
young population < 5 Myr (Ma dden et al.||2006| ). 



3.3. Supplemental data 

We supplement the Spitzer and LABOCA observations with 
2MASS J, H, K flux density estimates available on the 
NASA/IPAC Infrared Science Archive to describe the stellar 
contribution of the SEDs of the galaxies Mrk 1089, Haro 11 
and NGC 1705. We assume that these 2MASS flux densities 
are global measurements for the galaxies. For the UM 311 sys- 
tem, we perform the photometry and estimate the 2MASS flux 
densities from the images also obtained from the NASA/IPAC 
Infrared Science Archive. The IRAS broadband flux densities 
at 12, 25, 60 and 100 /mi are finally obtained through the 
NASA/IPAC Infrared Science Archive, t he IRAS Faint So urce 
catalogue, |b| > 10, Version 2.0 ( [Moshir & eTal) [1990), the 
IRAS catalogue of Point Sou rces, Version 2.0 (IPAC 1986) and 
Dultzin-Hacy an et al. ( |1990| ). They are added to supplement the 
existing data on dust emission. For Haro 11, the IRAS fluxes 
are consistent with the values of the IRS spectrum. For UM 311, 
IRAC data should probably encompass emission from both of 
the nearby galaxies and are thus not used in the analysis. 



3.2.3. IRS 

For the galaxy Haro 11, we supplemented our dataset with the 
mid-IR spectrum from the Spitzer Infrared Spectrograph (IRS) 
in order to better constrain the 5 /mi to 40 /mi range of the 
SED. We used public released data (Spitzer AOR key: 9007104) 
performed in staring mode. We used the Basic Calibrated Data 
products available in the Spitzer archive for the short-low resolu- 



4. Images 

4.1. Discussion on the morphology 

The 870 /mi LABOCA maps are presented in Fig. [2] The contours 
of the IRAC observations - 4.6 /mi for Haro 11 and 8 /mi for 
UM 311, Mrk 1089 and NGC 1705 - are superimposed on the 
LABOCA images. 
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Fig. 3. a) The positions of the short-high (green) and low-high (red) slits of the IRS observations superimposed on the Spitzer/3.6 
jim image of Haro 1 1 . b) The short-low and short-high spectrum of Haro 1 1 . 



NGC 1 705 - There is significant change in the galaxy mor- 
phology from NIR to FIR wavelengths. The bright 3.6 jim source 
is observed toward the location of the SSC while two strong 
MIR and FIR emission peaks appear at 5.8 jim and longer wave- 
lengths, offset from the SSC and coincident with the Ha maxi- 
mum (Fig. [T] and Fig. [2^). These two off-SSC emitting regions 
do not have bright stellar counterparts. The Spitzer 8 jim image 
shows the SSC in the middle of these two off-nuclear HII re- 
gions. The eastern region is the brightest source at 24 jim, with a 
flux density which is two times higher than the SSC flux density 
at this wavelength ( Can non et al. 2006 ). The three regions are 
also barely resolved in the 870 jim images. An offset peak is de- 
tected toward the west of the centre of the galaxy. Two faint 24 
jim sources are possible counterparts of the emission. 

Haro 11 - This galaxy is not resolved by Spitzer at wave- 
lengths greater than 4.5 jim. It is clearly detected but barely 
resolved with LABOCA (Fig. [2j)). The final images present ex- 
tended structures but this extension is not observed in the MIR 
images and was removed while calculating the submm flux to be 
conservative. 

Mrk 1089 - The 8 jim image (Fig. |2j:) clearly shows the 
interaction between the galaxies NGC 1741 (East) and Mrk 
1089 (West). A Spitzer color composition image of the complete 



Hickson group observed with Spitzer is described in Johnson 
et al.| ( |2007| ). The merging center of Mrk 1089 dominates the 



emission at MIR, FIR and 870 jim. Diffuse emission probably 
linked with intergalactic dust, presumably from galaxy interac- 
tions, is detected through out the whole region and was removed 
while performing the photometric measurements. 

UM 311 - The 4.5 jim image of the interacting field (Fig.|2jl) 
shows the spiral structure and arms of NGC 450 as well as the 
spiral companion UGC 807 and the three bright HII regions 
between these two spirals, one of which is the UM 311. We 
adopt the same numbering as |Moles et~aL] ( |1994| ) to describe 
these 3 HII sources and will call them "the UM 311 system". 
UM 311 is the region called 3 on the 4.5 jim image. This com- 
pact HII galaxy is the brightest source of emission at 24 jim. In 
the LABOCA image, the emission peaks toward the location of 
UM 311. The galaxy should dominate the SED of the UM 311 
system. 

4.2. Photometry 

We use the function aper of the library IDL astro of the NASA 
Goddard Space Flight Center to perform the aperture photome- 



try. For the Spitzer photometric data, the background is estimated 
using annuli just outside the boundaries of our galaxies. For the 
LABOCA images, the background significantly varies in inten- 
sity throughout the entire map. We calculated the flux densities 
in small circles in the immediate surroundings of the galaxy and 
averaged them to determine a local background for each source 
which was removed during the photometric calculation. 

To determine the flux densities at IRAC and MIPS bands of 
our sources and compare data between different wavelengths, 
the observations are convolved and regridded to a common res- 
olution. For NGC 1705 and Haro 11, the lowest spatial resolu- 
tion is the resolution of MIPS 160 jim (40"). We use convolu- 
tion kernels ( [Gordon et al.||2008] ) which convert a higher resolu- 
tion IRAC/MIPS point- spread function (PSF) to lower resolution 
IRAC/MIPS PSF using Fourier transforms. For the LABOCA 
data, we convolved the observations with a Gaussian kernel and 
regridded the images to the resolution of MIPS 160 jim. The 
photometric aperture chosen has a diameter of 144" for these 
two galaxies to encompass the entire emission of the galaxy. 

As Mrk 1089 was not observed at 160 jim, the four IRAC 
bands and the MIPS 24 and 70 jim observations were regridded 
to the resolution of LABOCA (18.2"). We choose an aperture of 
72" avoiding the flux arising from the other nearby companions. 
However, Mrk 1089 clearly dominates the whole system at MIR- 
FIR wavelengths. 

For UM 311, the resolution of the MIPS bands and the prox- 
imity of the 2 nearby HII regions makes it difficult to study 
UM 311 alone. We decide to perform the photometry of the 
UM311 system (encompassing sources 1, 2, 3 in Fig.[2]i). The 
2MASS, IRAC and MIPS 24 and 70 jim images were regrid- 
ded to the LABOCA resolution of 18.2" for the UM 311 system. 
In the 160 jim image, the broad PSF may cause emission from 
other sources outside the region around UM 3 1 1 to bleed into 
that region (Fig.[2jl). To be conservative, we choose a 54" aper- 
ture which encompasses the UM31 1 system and use the 160 jim 
observation as an upper limit when modelling the SED of this 
system. The procedure we use to model this system gives us so- 
lutions for the limits on the parameter space for UM 3 1 1 as if this 
interacting system were more distant causing these 3 sources to 
blend together. 

The IRAC measurements require aperture corrections to ac- 
count for the scattering of incident light in the focal plane arrays 
( [Reach et al.) [2005] ). These scaling factors of 91, 94, 71, 74% 
are respectively applied to the final integrated flux density esti- 
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NGC 1705 


Haro 1 1 


Mrk 1089 


UM 311 system 


2MASS 


1.25 jim flux density 


(in mJy) 


52 + 1.6 


13 + 0.3 


12.2 + 0.6 


12.8 + 1 


2MASS 


1.65 fim flux density 


(in mJy) 


50.9 + 2 


13 + 0.5 


13.9 + 0.9 


12.6 + 1 


2MASS 


2.17 jim flux density 


(in mJy) 


41.1 +2.4 


13.8 + 0.7 


12.2 + 1.2 


11.5 + 1 


IRAC 


3.6 jim flux density 


(in mJy) 


26.3 + 3 


22.5 + 2 


14.9 + 2 




IRAC 


4.5 jim flux density 


(in mJy) 


17.9 + 2 


28.9 + 3 


11.2 + 1 


5.5 + 0.5 


IRAC 


5.8 jim flux density 


(in mJy) 


16 + 2 


72.7 + 7 


33.3 + 3 




IRAC 


8 jim flux density 


(in mJy) 


17.9 + 2 


177 + 18 


92.5 + 9 


35.6 + 4 


IRAS 


12 jim flux density 


(in mJy) 




417 + 46 


106 + 21 




MIPS 


24 jim flux density 


(in mJy) 


52.7 + 5 


1930 + 193 


460 + 46 


151 + 15 


IRAS 


25 jim flux density 


(in mJy) 




2500 + 25 


579 + 29 




IRAS 


60 jim flux density 


(in mJy) 


868 + 61 


6880 + 41 


4010 + 40 




MIPS 


70 jim flux density 


(in mJy) 


1110 + 222 


4190 + 838 


3990 + 798 


1480 + 148 


IRAS 


100 jim flux density 


(in mJy) 


1610+100 


4730 + 640 


5430+ 181 




MIPS 


160 jim flux density 


(in mJy) 


1110 + 222 


1660 + 332 






LABOCA 


870 jim flux density 


(in mJy) 


114+17 


40 + 6 


67.4 + 10.1 


32.2 + 7 



Table 3. Integrated flux densities measured with 2MASS, IRAS, Spitzer/IRAC, Spitzer/MlPS and LABOCA. The 10% non-dust 
contamination was subtracted from the LABOCA flux densities. For Haro 1 1, no broadband data between 5 and 38 /im were used in 
the SED modelling since we add the IRS spectrum to constrain this range. 

* The available IRAS fluxes for UM 311 are not used as they incorporate the entire field which includes the two nearby spiral 
galaxies too. 

** The 160 jim resolution is too large to isolate the 160 jim flux of the UM 311 system alone. 



mates at 3.6, 4.5, 5.8 and 8 /im. Table [3] presents the IRAC flux 
densities corrected for the scaling factors. 

Finally, we compare the fluxes obtained from Spitzer obser- 
vations with the IRAS broadband fluxes. MIPS 24 jim and IRAS 
25 jim as well as IRAS 60 jim and MIPS 70 jim are consistent 
for the galaxies where both are available except for Haro 1 1 for 
which IRAS fluxes are systematically higher than MIPS fluxes 
due to the fact that IRAS fluxes were estimated using bigger 
aperture. The IRAS 25 jim flux is not used since we use the IRS 
spectrum as a constraint for that part of the SED. The IRAS 100 
Iim flux values seem to be high for Mrk 1089 and NGC 1705, 
perhaps encompassing more extended emission. We will see in 
§5 what this implies for the SED modelling. 



Sky Survey (843 MHz - |Mauch et al.] [2003). We derive the 
870 [im radio continuum contribution by extrapolation of the 
radio data tendancy (v -1 ). We find a radio contribution of 3% 
of the LABOCA flux for Haro 11. The radio observations in 
the literature are not sufficient to constrain the expected 870 
Iim radio contamination for NCG1705, UM 311 and Mrk 1089. 
We thus consider that the radio continuum contamination is 
of the same order for these other galaxies. Finally, we conser- 
vatively estimate the global non-dust contamination from po- 
tential CO(3-2) emission and radio continuum emission (syn- 
chrotron/bremsstrahlung) in the 870 jim band to be 10% for our 
four galaxies, similar to the submm contributions that Galliano 
et al. (2005), for example, determine for their dwarf submm ob- 
servations. 



4.3. Contamination of the 870 \im emission from non-dust 
origins 

Since we are interested here in the dust emission arising in the 
870 /mi band, we must consider and quantify possible contami- 
nation from non-dust sources and take these corrections into ac- 
count when modelling the dust SEDs. The CO(3-2) line can, in 
principle fall within the 870 jim band. While CO(1-0) has been 
a great challenge to detect in low metallicity galaxies (c.f. |Leroy| 
|et al.| |2005| ), we can not be sure that the smaller beam, higher 
excitation CO(3-2) could not be present here. CO observations 
were attempted in NGC 1705 without positive detections ( [Greve 
et al.||199"6| ). In the same way, CO seems to b e very faint in th e 



Hickson Group to which Mrk 1089 belongs (|Yun et aL]|1997|). 



Finally, for Haro 11, Berg vail et al.| ( |2000| ) found an upper limit 
for Lco(i-o) of 10 29 W. We can derive a n upper limit t o Lcq q-2) 
from the L C o(i-o) estimation using the |Meier et aL| ( |2001| ) re- 
lations for dwarf starburst galaxies: Lco(3-2)/Lco(i-0) is usually 
lower than 1. For the four galaxies, we conservatively estimate 
the CO (3-2) contribution to the 870 jim band to be ~ 5%. 

Additionally, we can expect contributions to the 870 /im 
fluxes from radio continuum emission (synchrotron emission 
and/or bremsstrahlung). Three radio fluxes were estimated for 
H aro 11 with the NRAO VLA Sky Survey (8.46 and 1.4 GHz 
- Condon et al. 1998) and the Sydney University Molonglo 



5. The SED modelling 

The SED is a complex tool to interprete. For a macroscopic re- 
gion of the ISM (HII region, cirrus), the SED summarises a wide 
range of physical conditions linked with the non-uniformity of 
its illumination and the variations of the dust composition. For 
an entire galaxy such as the dwarf galaxies we are studying here, 
the SED synthesises, on the contrary, the different components 
that constitute its global shape and emission: HII regions, molec- 
ular clouds, nebulae, diffuse ISM etc., that can not be studied 
independantly due to the lack of observational constraints. 

We want to construct the SEDs of our galaxy sample us- 
ing the fluxes presented in Table [3] and the IRS spectrum for 
the galaxy Haro 1 1 as constraints in order to quantify elemen- 
tary quantities of these galaxies such as their mass of PAHs or 
their total mass of dust. The stellar contribution to the SED is 
constrained by the 2MASS bands and IRAC 3.6 jim. The NIR 
to submm wavelengths of the SEDs, signatures of the dust and 
physical conditions of the galaxies, are constrained by IRAS, 
Spitzer and our new LABOCA data (Table [3]). No radio com- 
ponent is taken into account as a model constraint but the radio 
contribution is removed from the 870 jim flux (see § 4.3). 

The model we use to fit the SEDs is a simplified version 
(due to the smaller nu mber of observational constraints) of the 
Galliano ^et al.| ([2008) model. We would first like to note that a 
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SED depends on the intensity and the hardness of the Interstellar 
Radiation Field (ISRF) that the ISM experiences, on the mass 
fraction of the dust species (silicate and carbon grains, PAHs 
etc.) and on the distribution of grain sizes. Thus, modelling a 
SED requires making a priori assumptions on the ISRF and on 
the global properties of the dust, assumptions that we will de- 
scribe in the following paragraphs. 

We assume that the source of excitation of the dust is the 
ISRF which we choose to have the spectral shape of the ISRF 
found in the Galactic diffuse ISM ( [Mathis et aL][T983] ), although 
we test the influence of this spectral shape by trying a variety of 
forms (see § 6.4). The ISRF intensity will be scaled using a fac- 
tor U (defined in |Draine & Li||2007] ), with U=l corresponding 
to a normalization to the local solar neighbourhood value of 2.2 
x 10" 5 Wm" 2 . 

We assume that the sources of IR emission are dust and 
old stars. We suppose that the dust composition is homoge- 
neous throughout the galaxy. For silicates, graphites and PAHs, 
we ado pt the co mposition and size d istribution of |Zubko et al. 



P004| ) (Table [4|). [Zubko et al.| ( |2004| ) assume that the dust par- 
ticules are PAHs, graphite and silicate grains and that the ISM 
has a solar abundance. The assumed optical properties of these 



grains are taken from Draine & Li ( 2 007 1 PAHs), |Laor & Draine 
(1993, graphites) and Weingart ner & Draine ( 2001| silicates). 
The total mass of dust (M^) represents the first parameter of 
our model while the PAH component requires the introduction 
of two other parameters, the ionised PAHs-to-neutral PAH ratio 
(fpAH+) and the PAH-to-total dust mass ratio (Xpah) normalised 
to the Galactic value of 0.046, to be described. 

We adopt the prescription of |Dale et al.| ( |2001| ) to relate the 
dust mass exposed to a given intensity dMd ust (U) to the different 
heating environment intensities (U) to which dust is exposed. 

dM dust (U) oc U~ a dU, U min <U< U max (1) 

This prescription is flexible enough to describe dense and 
diffuse media. The simple power law relation leads to the def- 
inition of the index a, which represents the contribution of the 
different local SEDs exposed to a certain intensity U. We intro- 
duce the three quantities a, the minimum heating intensity U m ;„ 
and the maximum heating intensity U max as free parameters in 
our modelling. 

Since the stellar component contributes to the NIR part of 
the SED, we finally add it from the dust SED by fitting a stellar 
spectrum. This spectrum is synth esised using the st ellar evolu- 
tion code PEGASE (Fio c & Rocca-Volm erange , 1997 ) assuming 
a Salpeter Initial Mass Function (IMF). The stellar population 
is considered to have undergone an instantaneous burst (5 Gyr 
ago) and the initial metallicity is assumed to be solar (Z=Z©). 
The mass of stars (Goldstar) is introduced as a parameter of our 
modelling. 

In summary, the free parameters used in the modelling are: 



and near-IR data are first used to constrain the stellar radiation 
of the model. We then compute the temperature distribution of 
the dust grains heate d by the absorption of this ste llar radiation 
using the method of Guhathakurta & Draine (1989). Synthesised 
spectra are finally computed for each of the silicates, graphites, 
neutral and ionized PAHs, considering that grains are stochasti- 
cally heated. The sum of these discrete contributions leads to a 
global SED model of the galaxy. 

The interactive fitting stops when deviations from the obser- 
vational constraints are minimized. The x 2 minimisation algo- 
rithm is based on the Levenberg-Marquardt methods. To prevent 
X 1 from being influenced by the density of points or dominated 
by the highly sampled MIR spectrum in the case of Haro 1 1 , we 
weight each data point depending on the density of points around 
its wavelength. A L° bs (Ai) represents the error on the luminosity 
at a given wavelength. 



X 



z(- 



2At 



L° v bs (Ai)-L v Ui) 



ALy(Ai) 



(2) 



We will call this first model the "fiducial" model. 

Our model presents some limits mainly linked to the assump- 
tions we made to simplify the mo del. Indeed, we choo se the ra- 
diation field shape of the Galaxy (Mathis et al. 1983] ), a profile 
which could be different in low-metallicity environments. The 
choice of a radiation field profile usually affects the PAH and 
small grains mass derived from the modelling. The influence 
of this parameter is disc ussed in §6 . 4. Ano ther strong assump- 
tion was in the use of the Dale et al. ( |2001| ) prescription (Eq.[T]). 
This formula directly relates the cold and the hot regions of the 
galaxy, an assumption which is not valid if cold dust is physi- 
cally residing in different regions than star forming regions. The 
fact that a separate cold dust component is required for so me of 
our galaxies (see § 6.2) may be a sign of the limit of the Dale 
|et al.| ( |2001| ) prescription. Nevertheless, our simple model has 
the sufficient level of complexity (e.g. detailed dust properties) 
to accurately derive the global properties we want to study, es- 
pecially as regards to our small number of observational con- 
straints for each galaxy. Adding other components such as dif- 
ferent phases, clumps or more complex geometries would have 
lead to an overinterpretation of the observations. 

An international conference was held (2005) to discuss 
other SED models. We refer the readers to the reviews of this 
conference in|Popescu & Tuffs (2005 ). 



6. Analysis 

6.1. Results of the fiducial SED modelling 



M dust 

/PAH 
fpAH+ 

a 

U min 
f^ max 
Goldstar 



total mass of dust 

PAHs-to-dust mass ratio (normalised to the 
value) 

ionised PAHs-to-total PAHs mass ratio 



The SED models are presented in Fig.[4j along with the observa- 
tions. Haro 1 1 observations are displayed corrected to rest wave- 
G a i ac ti c lengths taking into account its redshift. The Haro 1 1 IRS spec- 
trum, from 5.1 yum to 37.2 /mi, provides additional constraints 
to better describe the PAH properties as well as the slope of the 



index describing the fraction of dust exposed to a given risin g MIR continuum. We perform SED models using the 870 
intensity " ^ m observational constraint, and then compare to the resulting 

minimum heating intensity parameters when we remove the 870 fim data to study the in- 

maximum heating intensity fluence of the submm data on the parameters of the modelling, 

mass of old stars especially the dust mass of the galaxy (Table [5]). To quantify the 

errors on our dust mass determinations, we produced a grid of 
500 randomly modified observational constraints, allowing the 
observational fluxes to vary within their error bars, following a 



The modelling is an iterative process in which we assume 
an initial dust grain distribution (size and composition). Optical 
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PAHs Graphites Silicates 



Minimum Size (//m) 
Maximum Size (pirn) 
Mass density (g cm -3 ) 



3.5 xlO" 4 3.5 xl(T 4 3.5 xlO" 4 
5 xlO" 3 0.33 0.37 
2.24 2.24 3.5 



Table 4. Size range (Zubko et al. 2004 ) and mass densities of the three dust grain components. 



Mrk 1089 



_l 





10 100 
Wavelength (yum) 



1000 



10 100 
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Fig. 4. SED models of Mrk 1089, the UM 311 system, NGC 1705 and Haro 11 using the fiducial model. The SEDs are plotted in 
black. Observational constraints (listed in Table[3]) are superimposed (filled circles). The green and red lines respectively distinguish 
the stellar and the dust contributions. The dashed black lines present the SED models of our galaxies obtained when the LABOCA 
constraint is not used in the modelling. The open circles represent the expected modeled fluxes integrated over the instrumental 
bands. When the error bars are not shown, the errors are smaller than symbols. Note that the IRS MIR spectrum used in the 
modelling is overlaid in orange for Haro 11. For the UM 311 system of 3 compact sources, the 160 yum flux is an upper limit since 
it was calculated with a 40" aperture. The different SEDs represent the possible SED models that fit the observational constraints 
with good accuracy. 



Gaussian distribution around their reference value. Before enter- 
ing into the details of each galaxy, we can already notice that the 
fpAH, which is normalized to the Galactic value, is always infe- 
rior to 1, which is consistent with the fact that low-metallicity 
galaxi es usually show a PAH deficit compared to du s tier galax- 
ies (|Madden| [20051 [Engelbracht et al.| [2005) [20081 |Wu et al.| 
2007[). For NGC 1705, 11.3 /mi PAH emission was neverther- 



less detected with IRS spectroscopy in the most luminous of the 
two dust emission peaks (Dl) but not in the second off-nuclear 
region nor in the SSC ( [Cannon et al. 2006 ), so most of the 
PAHs should reside in Dl. Finally, a is quite similar for the 



three galaxies, with Haro 11 having the smaller a. Dale et al. 
( 200 51 [2007]) found that or was typically 2.4 for SINGS galaxies 
while |Draine et al.| ( [2007 1 ) fixed it at 2 in their SED modelling. 
Our values are smaller than 2 for NGC 1705 and Haro 11. Note 
that for a given U max , galaxies with smaller a usually contain 
more strongly heated regions and therefore are more likely to be 
associated with intense star formation activity. 

Mrk 1089 - In spite of the lack of 160 fim observations, we 
are able to fit the observations accurately both with and without 
the LABOCA observations. Notice how flat the peak of the SED 
is in the FIR (Fig. [4]). The slope of SINGS galaxies is positive 
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Table 5. Model parameter results derived from our fiducial SED modelling of Haro 11, Mrk 1089 and NGC 1705 and dust-to-gas 
mass ratios implied by the dust mass. Md ust represents the total mass of dust of the galaxy in solar mass, f P AH the PAH-to-total dust 
mass ratio, a the power law index describing the illuminating radiation fields, U min and U max the minimum and maximum heating 
intensities. Finally, Mm represents the HI mass and Mdustf^m the dust-to-HI mass ratio. Bold numbers are the dust mass and the 
dust-to-gas mass ratio values of our prefered SED for the galaxy Mrk 1089. For Haro 1 1 and NGC 1705 where the model does not 
fit correctly our observed constraints, a cold 10 K dust component is added (see § 6.2). 

* For NGC 1705, we use the HI mass associated with the aperture chosen for the dust mass determination to calculate the Mdust f^m- 
The HI mass reported here is thus 80% of the the total HI mass (5.1 x 10 7 M o ) of the entire extension of HI. 
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Table 6. Model parameter results derived from our fiducial SED modelling for the UM3 1 1 system. We perform the modelling using 
the upper limit of the 160 yum flux and decrease it in steps of 10% down to 50% to study the evolution of the model parameters. 
These parameters are: Md us t the total mass of dust of the galaxy, f P AH the PAH-to-total dust mass ratio, a the power law index 
describing the illuminating radiation fields, U min and \J max the minimum and maximum heating intensities. 



between 60 and 100 yum whereas the slopes of the dwarf galaxies 
here (Mrk 1089 but also in NGC 1705 and Haro 11) are flat or 
descending in that wavelength range. This rather flat form of the 
SED is not characteristic of our Galaxy nor is it common in other 
more metal rich galaxies, but is often seen in the SEDs of active 
low metallicity galaxies, such as IIZw40 and He2-10 ( Galliano 
20031 12005] ), in the Larg e Magellanic Clou d ( |Bernard| 
20081 or in the SEDs of |Dale et al] ( [2007] ). There are 



et al. 




et al. 





two possible explanations for this: the lower metal abundance 
and the decrease in dust attenuation result in a higher overall in- 
terstellar radiation field and/or the increased abundance of very 
small stochastically-heated grains can inflate the continuum on 
the Wien side of the SED. We find a total dust mass of 5.12 x 
10 7 Mq (± 20%). This mass is a factor of 1.7 larger than that 
we would calculate without submm constraints. This galaxy is 
the least affected by the submm constraint. Nevertherless, this 
increase is superior to the error bars estimated for this galaxy. 

The UM 311 system - The 160 yum beam is too large to isolate 
the UM 311 system alone and, without this data, we lack suffi- 
cient observational constraints at FIR wavelengths to accurately 
determine the shape of the SED peak for the system. Instead 
we measure the 160 yum flux in a bigger aperture (80") than the 
one chosen to measure the fluxes in the other bands (54") and 
use the 160 yum flux as an upper limit. To explore the range of 
parameter space constrained by the 160 yum upper limit, we de- 



crease the 160 yum flux in steps of 10% down to 50% and find 
the range of shapes with associated parameters that could fit our 
observational constraints (see Fig. [4] and Table [6]). We also pro- 
duced a grid of 200 randomly modified observational constraints 
to study the spread of the SED models induced by these uncer- 
tainties. These tests enable us to establish a range of possible 
dust masses for this galaxy system. We conclude that the dust 
mass of the system should reside between 5.2 x 10 6 and 1.1 x 
10 7 M Q , and a should reside between 1.79 and 2.13 (Table [6}. 
The various SED models obtained include both solutions that 
maximise the contribution from the hot dust with a FIR peak at 
about 50 jim and solutions that, on the contrary, favor a large 
dust mass, peaking at longer wavelengths. 

We convolve the 8 yum IRAC images of the UM 311 system 
to the resolution of MIPS 24 /im and calculate the luminosities 
of the 3 different substructures at these wavelengths. The lumi- 
nosities in y L y of region 1 and 2 are respectively 1.5 and 2.9 
times lower than that of UM 311 (region 3) at 8 /mi and respec- 
tively 2.6 and 8.4 times lower at 24 yum. Thus, we obtain dif- 
ferent 8/24 luminosity ratios for the 3 different regions, respec- 
tively 0.36 and 0.61 for region 1 and 2, and 0.21 for UM 311 
which compares with the value of 0.24 obtained for the whole 
UM 311 system. These ratios imply that region 2, which has the 
lowest 24 yum flux, may be relatively less active than the others, 
if we use the 24 yum flux as a good indicator of the star formation 
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activity of the region (|Helou et alT] |2004[ [Relano et al.[ 2007) 



On the contrary, UM 311, having the highest 24 jum continuum 
should present the strongest star formation activity of the sys- 
tem. Moreover, even if the system is not clearly resolved at 70 
/mi, its emission clearly peaks at the location of UM 311. This 
leads to the conclusion that the galaxy UM 311 should dominate 
the system at MIR and FIR wavelengths and the global dust SED 
of the UM 311 system. The galaxy should account for the major 
part of the hot dust mass of the system. If more distant, the 3 
objects would blend, like the 3 distinct nuclei of Haro 1 1 which 
can be resolved by HST but not resolved with Spitzer. 

NGC 1705 - Fig. [4] presents the fit which gives the lowest 
X 2 for this galaxy. The model fails to fit the IRAS 100 yum data 
when we introduce the LABOCA constraint but does fit this con- 
straint when the SED is modeled with constraints up to 160 jim. 
We test the influence of the IRAS 100 jim flux on the fit per- 
forming SED models with and without this data point but with 
the submm constraint. The global shape of the SED does not 
change and t he dust masses derive d in these 2 cases vary by 
less than 5%. |Cannon et al.| ( |2006|) estimated th e dust mass of 
NGC 1705 using the models oT |Dale & Helou] p002) and |Li 
& Draine| ( |200l) , obtaining respectively ( 3.8 + 1.9) x 10 5 Af ? 
and 7 x 10 4 M Q (± 50%). The model of ]Li & Draine] ([2001 1 



combines laboratory studies and astronomical observations and 
fits the MIPS fluxes very well. The dust mass deduced from the 



Li & Drame| ( |2001| ) model, accounting for its 50% uncertainty 



and the different distance they used (5.1 Mpc), compares well 
with our results obtained without the submm constraint: 3 x 10 4 
M Q (Table[5|. Including the submm LABOCA constraint, we ob- 
tain a dust mass of 1.3 x 10 6 M , which is 50 times higher than 
our results without the submm data. In any event, our fiducial 
model, even using the 870 jim observations, is not satisfactory 
and gives a poor^ 2 value. For this galaxy, the uncertainty in the 
mass estimate obtained using the submm constraint is mostly 
due to the flattening FIR peak and to the elevated 870 jim emis- 
sion, for which the model has difficulty to find a consistent so- 
lution. Considering that the error in the dust mass reaches more 
than 50%, we consider this first dust mass estimate obtained with 
our fiducial model uncertain and not satisfactory for this galaxy. 

Haro 11 - The SED of Haro 1 1 is very striking in that it peaks 
at very short wavelengths - 36 yum, highlighting the extreme na- 
ture of the young starburst. This is also evidenced by the high 
global value of the Ne[III]/[NeII] ratio (r atio > 1), already seen 
in other low metallicity dwarf galaxies (Madd en" et al.| |2006| ), 
which is normally an indication of the hard interstellar radia- 
tion field dominated by a young (< 5 Myr) stellar population. 
The peak of the IR SED compares to that of the 1/4OZ galaxy 
SBS0335-052 flHouck et al.| [2004] ) or IC10 NW (Parkin et al. 
2009 in prep). The IRS spectroscopy provides a tight constraint 
on the slope of the hot MIR dust continuum emission as well as 
for the fpAH+- The model does not clearly fit all the details of the 
IRS spectrum due to the lack of sophistication and flexibility in 
the dust properties. More complex modelling would be required 
to perfectly model the different features of the spectrum but this 
does no t affect the following con clusions on the dust mass pa- 
rameter. Engelbra cht et al.|(2008| ) previously estimated the total 
dust mass of the galaxy to be 6.2 x 10 6 M Q using Spitzer ob- 
servations and a distance of 87 Mpc. This value is consistent 
with the dust mass we obtain without submm constraints for a 
distance of 92 Mpc: 7.2 x 10 6 M . Note that when the IRS spec- 
trum is not used for the modelling but only IRAC broabands, the 
dust mass is estimated to be 2 x 10 8 M© with a poor-fit (J 2 = 
24). Thus the IRS spectrum greatly influences and constrains the 
global shape around the peak of the SED. 



The large uncertainty in the dust mass estimated from the 
modeled SED (the error reaches 40 % for this galaxy) is due 
to the excess emission at submm wavelengths that the model 
does not fit (Fig. [4]). Such submm excess was already observed in 
other galaxies using submm instruments such as SCUBA on the 
James Clerk Maxwell Telescope (JCMT) ([Galliano et al.| [2003] 
[2005l|Bendo et al.| [20061 |Marleau et al.||2006|) and usu ally leads 
to different interpretations. |Galliano et al. (2003 , 2005 ) proposed 
that the submillimeter excess in the dwarf galaxies that they ob- 
served could originate from 10 K dust. We decide to explore this 
hypotesis. Our fiducial model does not seem to be sufficient to 
explain the fluxes we obtained at submm wavelengths for the 
two galaxies NGC 1705 and Haro 11. Two separate thermal dust 
models may be more successful to describe the 60-160 fim and 
the 870 yum emission separately for these two galaxies. 

6.2. Refinement: adding a cold dust component 

In an attempt to improve the model results, we refine our fidu- 
cial SED model for the two galaxies Haro 11 and NGC 1705 by 
adding a thermal cold dust component with a flux, S y , charac- 
terised by a modified blackbody: 



S v oc^B v (T d ) 



(3) 



where B v is the Planck function, T^ is the dust temperature and 
the emissivity coefficient ft. As we do not have information on 
the properties of these cold grains, we choose to describe them 
by extrapolating the properties of graphites grains and assume 
an emissivity ft of 1 or 2. We adopted the values p=2.5 x 10 3 kg 
m~ 3 , /io=100 jim and Qo/ao=150 cm -1 for the mass density , the 
reference wavelength and the absorption coeffcient at Aq of the 
grains respectively. We note that the absorption opacity chosen 
may slightly deviate away from the |Li & Draine| ([2001) value 
when the emissivity is 1 since they usually assume that j3=2. We 
nevertheless assume that their absorption opacity will still be 
applicable for j3= 1 . 

Indeed, there is a bend in the SED near 450 fim. Due to the 
physical parameters of the grains behaving like a modified black 
body in the submm, having a natural bend on the Rayleigh- Jeans 
slope seems complicated unless you suddenly change the emis- 
sivity of the grains at 450 /mi. In § 6.2, we discuss other possible 
explanations beyond the cold dust hypothesis. 

Since the supplemental black body component is only con- 
strained by the 870 yum flux, we cannot solve for the temperature 
and ft independently. For this exercise, we fix the temperature to 
10K and fit the SED using a /3 = 1 or 2 to test the hypothesis of a 
submm exce ss of cold dust ( Galliano 'et al. [ |2003] |2005 [ |Mar leau 
|et al.||2006| ). The new SED produced is obtained by adding the 
fiducial SED model and the cold blackbody. Chi- square values 
are estimated from the deviations of the new SED model to the 
observational constraints. In this revised model, the warm dust 
component is described by the observational constraints cover- 
ing the FIR wavelengths to 160 /mi while the cold component is 
constrained by the 870 jum observations. 

The dust masses derived from each SED model for Haro 1 1 
and NGC 1705 using J3=l and 2, as well as the dust masses de- 
rived from the fiducial models, are given in Table [7] Including a 
10 K dust component, we find dust masses constrained by j3 = 
2 to be 6 and 8 times greater than those with J3 = 1 (Table [7]) for 
Haro 11 and NGC 1705 respectively. Indeed, since a 1 OK black- 
body modified with a A~ l emissivity is flatter on the millimetre 
tail, the model fits the data with less mass than when a steeper 
emissivity coefficient of 2 is used. For Haro 11, the best-fit model 
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Fig. 5. SED models of NGC 1705 and Haro 11, adding a cold dust component of 10K with an emissivity coefficient J3 of 1 and 2 to 
our fiducial model. Observational constraints are superimposed with filled circles. The IRS MIR spectrum is overlaid in orange for 
Haro 11. The green, red and blue lines represent the contributions of stars, warm and cold dust respectively. The red lines is close 
to the model generated when submm is omitted in the SED fit. The new models fit well the FIR peak and the 870 fim flux. 
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Table 7. Model parameter results for the SED model introducing a cold dust component at 10K and an emissivity J3 of 1 and 2 to our 
fiducial model for the galaxies Haro 1 1 and NGC 1705. Results are compared to the ones of the fiducial modelling. Md ust represents 
the total mass of dust of the galaxy. This mass includes the cold (10K) dust mass for Haro 11 and NGC 1705. Md ust /M H [ represents 
the dust-to-HI mass ratio and fpAH the PAH-to-total dust mass ratio. Bold numbers are the dust mass and D/G values of our prefered 
SED models. 



prefers 0=1, giving a total mass of dust (warm and cold) of 2 x 
10 7 M Q , although the success in using 0=1 or 2 is not very dif- 
ferent in terms of the reduced chi- square (J 2 ). For both galaxies, 
the^ 2 values for the models including a cold component indicate 
better fits than the fiducial models. The SED models including a 
cold dust component are presented in Fig [5] Comparison by eye 
with the models presented in Fig. [4] show obvious improvements 



in the model fits to the observations. The value of significantly 
influences the dust masses derived from the SEDs and thus the 
dust-to-gas mass ratio (D/G) estimate of the galaxies (see section 
§5.3). We also note that the abundance of PAHs relative to the 
total dust mass of these galaxies is much lower than the Galactic 
value of 4.6 x 10" 2 : a factor of 15 lower for Haro 1 1 and a factor 
of 50 times lower for NGC 1705 (Table 0. 



M. Galametz et al: The SED of dwarf galaxies using LABOCA 



15 



For Haro 11, increasing the temperature of the cold dust with 
J3=l does not significantly influence the x 2 value until the tem- 
perature of the cold dust reaches ~20K. For NGC 1705, when 
we try to increase the temperature of the cold dust component, 
the 160 jim observations becomes difficult to fit and x 2 values 
increase. Thus, at least for NGC 1705, the very cold dust 10K 
dust provides a better fit, while for Haro 11, somewhat higher 
temperatures could also work. However, we need more observa- 
tional constraints to pin down the precise temperature of the cold 
dust component. 

As mentionned previously, Galliano et al. ( 2003 , 2005) pro- 
posed that the submillimeter excess in the dwarf galaxies that 
they observed could originate from 10 K dust with a dust emis- 
sivity index j3=l and found that dust would account for 40-80% 
of the total dust mass in each of these galaxies. We choose to 
follow the same assumption b ut other explanations can be pro- 
posed. Lisen feld et al.| ( |2001| ) proposed that the submillimeter 
excess would originate from hot (~100 K) dust with a dust emis- 
sivity i ndex J3= 1 and the temperature fluctuations of very small 
grains. |BencTo et al.| ( [2006] ) studied NGC 4631 at 450 and 850 
//in with SCUBA and found that the 850-1230 fim emission ex- 
ceeds what would be expected from thermal emission but the 
scenario of a very cold dust component which might explain 
this excess was rejected due to the high dust-to-gas mass ratio 
it would imply. In fact, submm excess is not always detected 
in low-metallicity galaxies. Parkin et al. 2009 (in prep) investi- 
gate the very well sampled SED of the low-metallicity galaxy of 
the Local Group IC10, using Spitzer and ISOCAM observations 
combined with JCMT/SCUBA observations. They model the 
SED of the 2 main SF regions only and found no submm excess. 
The observations, however, do not cover the entire gala xy, pre- 
venting analysis on the global scale. Draine et al. ( 2007| ) present 
SED models of a large sample of 65 galaxies. For 17 galaxies for 
which they have SCUBA data at 850 /mi, they fit SED models 
with and without this submm constraint. Five of their 17 galaxies 
show an increase in the dust mass when calculated using SCUBA 
constraints whereas 5 other galaxies show a decrease of this 
dust mass. In fact, their dust masses obtained with and without 
SCUBA data agree to within a factor of 1.5 for 11/17 cases and 
to within a factor of 2.2 for all cases. They concluded that their 
dust models do not require cold (<10K) dust to account for their 
submm fluxes. However, they caution that some of the SCUBA 
galaxies are not completely mapped or are not taken in scan map 
mode. The difficult data processing might have oversubstracted 
diffuse, extended emission for these galaxies. Thus assumptions 
they make for complet submm flux values may require caution 
in using these data. Their sample studied with submm observa- 
tions contains only metal-rich galaxies that may show different 
dust properties and dust temperature distributions than the low- 
metallicity galaxies we present here. Further studies would be 
necessary to study the presence of submm excess and the influ- 
ence of submm constraints on the SED model with respect to the 
metallicity of the galaxy. 

Instead of a cold dust component fitting the submm excess, 
we also test the formalism presented in Drain e et al. (2007 ) and 
try to fit the excess with an extra ISM component heated at 
U=U m ;„ for our two galaxies. This addition, in both cases, leads 
to a very cold component with U m ;„<10~ 2 due to the fact that 
the extra component tends to fit the submm excess. The masses 
of dust derived from the SED models in this case are 3.3 x 10 6 
M Q for NGC 1705 and 1.2 x 10 8 M Q for Haro 11, thus 19 and 6 
times higher, respectively, than the mass of dust derived from a 
SED model using a cold blackbody of 10K with fi=l, unrealistic 
values due to the D/G mass ratios they imply. 



Additionally, a j3= 1 .5 , which is closer to amorphous carbons 
(Rouleau & Martin, 1991), might give a better fit and reduce the 
cold dust mass but it would also imply that the FIR-mm dust 
mass was dominated by carbons and not silicates , which is con- 
trary to our usual relative dust mass concepts. 

Other assumptions can be invoked to account for the submm 
excess. For example, underestimating the contributions to the 
870 /mi flux values could have important consequences on the 
values of the deduced dust mass. Direct measurements of the 
CO(3-2) line and more mm radio observations would place 
greater confidence in the dust mass determinations. The submm 
excess could also be explained as a change in dust emissivity 
in lieu of the cold dust hypothesis. Models of |Meny et al.| ( |2007] ) 
modify the dust optical properties to find an effective decrease in 
the submm emissivity index as the dust temperature increases as 
suggested by the o bservations of |Dupac et al.| ( [2003] ). However, 
|Shetty et"aL] ( |2009| ) express caution in the inverse temperature 
- p interpretation, showing that flux uncertainties, especially in 
the Rayleigh- Jeans regime, can affect the results for the SED fits 
as far as temperature and emissivity are concerned. Although we 
do not reject the possibility of a modification of the properties of 
dust with temperature, here we explore the cold dust hypothesis, 
since it enables us to investigate its consequences on the global 
properties of the galaxies. 

6.3. Dust-to-gas mass ratios 

We calculate the D/G of our galaxies and compa re them to es- 
timate s from the chemical evolution model of Galli ano et al.l 
(2008). Values are presented in Table[5]for Mrk 1089, in Table[7] 
for Haro 1 1 and NGC 1705. Since our SED for UM 3 1 1 contains 
3 individual targets, we do not consider it here in the discussion 
of the D/G. 

NGC 1705 - This galaxy is surrounded by a very large HI 
envelope up to 10 times larger than the optical extension (see 
Fig. [I]), which is more often the case for dwarf galaxies in com- 
parison with spiral galaxies. From the HI image of NGC 1705, 
we estimate that only 80% of the total HI flux is contained in 
the aperture we choose to derive our dust mass. We then con- 
sider that only 80% of the HI mass, i.e 4.1 x 10 7 M Q , should 
be used for the D/G estimate. The D/G of NGC 1705 estimated 
with submm observations is ~ 4.1 x 10" 3 with our required very 
cold dust component using a ft emissivity factor of 1 (Table [7]). 
This value compares to that expec ted by the chemical evolution 
model used in |Galliano et alJfcOOS) (2 x 10" 3 < D/G < 5 x 10" 3 
). The solution for J3=2, obtained with a slightly higher^ 2 , leads 
to a D/G of 2.8 x 10" 2 , which is a factor of 8 higher than that 
of J3= 1 . This high D/G value would be difficult to reconcile with 
chemical models. We prefer the model with a cold component 
using £=1 for NGC 1705. 

Haw 11 - The D/G estimated with our SED model is ~ 0.2 
when a very cold dust component of J3=l is used (Table [7]). If 
we were to use an emissivity index of 2 for the very cold grains, 
the D/G derived is even higher and not physical considering the 
lack of metals in the galaxy. For this reason, in spite of the sim- 
ilar^ 2 values for the two different fits, we prefer the fit for f3=l 
for the cold dust component. The HI mas s used for our calcula - 
tion is an upper limit of 1 8 M© given in |Bergvall et al.|p000| . 
|Bergvall & Qstlin| ( 2002 ) already note from their VLA obser- 
vations that Haro 1 1 is dramatically deficient in neutral hydro- 
gen. For such a low metallicity (12+log(0/H) = 7.9), a D/G of 
~ 10" 3 should be more likely considering chemical evolution 
models. Our high values for the D/G, constrained by the new 
submm data, are not easily explained with current models. For 
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the same galaxy, Gallian oet al.| p008) estimated the D/G -1.6 
x 10" 2 (without submm constraints) but this ratio was qualified 
as uncertain due to the poor HI measurement for this galaxy 
which could be linked to the difficulty of accurately measuring 
HI at such distances. Indeed, in Gordon & Gottesman ( 1981|), 
the mean HI mass for BCGs of the size of Haro 1 1 should be ~ 
8.1 x 10 8 M Q . If we consider that the existing HI observations 
provide a reliable upper limit on the mass of atomic hydrogen 
( [Bergvall et al.| |2000| ), the large estimated dust mass may sug- 
gest that a significant fraction of the total gas mass is in forms 
other than atomic. A major fraction of its gas should be in molec- 
ular or ionised form. An upper limit of only ~ 10 9 M was es- 
timated from CO observations for the mass of molecular gas in 
this galaxy ( Bergvall et al. 2000| ). The total amount of observed 
molecular + atomic gas (all upper limits) still can not lead to a 
D/G of 10" 3 for Haro 11. In fact, CO can be a poor tracer of 
molecular gas in low metallicity galaxies due to the high excita- 
tion and the density of these environments and the small filling 
factor of the molecular clouds. Self- shielding can be extremely 
efficient for H2 in regions where CO is photodissociated. Thus, 
there could be a n on negligeable amount of molecular gas no t 

"fT995l|Madden et al " 



traced by CO (e.g Poglitsch et al. 



1997). 



In order to reach the D/G value expected by chemical models, 
a large molecular gas mass on the order of 10 10 M© would be 
required to account for the "missing" gas mass, that is to say an 
order of magnitude more molecular gas than that deduced from 
current CO upper limits. Molecular hydrogen co uld be embed- 
ded in [CII] emiting envelopes as described in Madden et al. 
( |1997| ) for the low-metallicity g alaxy IC10. Th e high L C ///L C o 
value (4 x 10 5 ) deduced by |Bergvall et al.|(|2 000) and also seen in 
other low metallicity galaxies (Madden, 2000) is coherent with 
this theory. 

Mrk 1089 - Our SED model estimate of the D/G leads to 
a value of 1.9 x 10" 3 . For a galaxy with the metallicity of 
Mrk 10 89 (~ 12+log(Q/H)= 8 .0), the chemical evolution model 
used in Gallian oet al.| ([2008 ) which links the D/G with metal- 
licity, leads to an estimated D/G of ~ 10" 3 , consistent with the 
value we obtain to within a factor of 2. 

To summarize, the SED results we will discuss in the fol- 
lowing sections are those obtained with our fiducial model for 
Mrk 1089 (i.e., no additional modified black body) and the SED 
models which include a cold dust component of 10K with J3=l 
forNGC 1705 and Haro 11. 



6.4. Robustness of the results with the assumed radiation 
field 

The modelling scheme already presented in this paper, does not 
constrain the form of the global ISRF of the galaxies and the 
results should be dependent on the shape of the average radi- 
ation field to which the dust grains are exposed. The SED re- 
sults p resented here are o btained using the Galactic shape of the 
ISRF (Mathi s~et al.| |1983| ), but this may not be accurate in our 
case for metal-poor galaxies. Modeled ISRFs of dwarf galaxies 
show that dwarf galaxies can certainly have h arder globa l av- 
erage ISRFs than that of the Galaxy ( [Madden et"aL| |2006| ) - an 
effect attributed to the lower dust attenuation in the low metallic- 
ity ISM and, as a consequence, the larger mean free path length 
of the ionising photons. As we are modelling galaxies on the 
global scale here, we test several ISRFs which have been deter- 
mined for low metallicity galaxies to quantify this dependence. 
We use synthesised radiation fields of four dwarf galaxies as a 
reference: He2-10, IIZw40, NGC1140 and NGC1569. The pro- 
cess to synthesized the ISRFs uses the stellar evolution synthe- 



sis code PEGASE ( |Fioc & Rocca-Vol merange 1997| ) and the 
photoionisation model CLOUDY ( [Ferland 1996| ). We refer to 
Galliano et al. ( 2003 , 2005) for more explanations of this process 
and a description of the synthesized ISRFs of He2-10, IIZw40, 
NGC1 140 and NGC1569. We furthermore produce a grid of dust 
models for the ISRF of a young and non-ionising cluster, created 
just after an instantaneous burst, with a Salpeter IMF. Even if the 
shapes of the resulting SEDs for our galaxies are somewhat mod- 
ified using different forms of harder radiation fields, we find that 
the masses of dust derived from these different models and those 
determined using the Galactic ISRF, differ by less than 10%, 
which is within the error bars of the dust masses we estimate 
(Table[8]). Indeed, the shape of the radiation field essentially con- 
trols the emission of out-of-equilibrium grains. Increasing the 
hardness of the radiation field increases the maximum tempera- 
ture that small grains reach when they fluctutate. Consequently, 
it produces a short wavelength excess of the grain spectrum as 
reflected in the resulting modified SEDs. However, this excess is 
compensated in the global model of § 5, by lowering the weight 
of the high intensity regions. Since these hot regions do not con- 
tribute significantly to the total dust mass, the latter does not 
strongly depend on the shape of the ISRF. 

6.5. Distribution of the dust temperatures 

To get an idea of how the dust mass is distributed as a function 
of the temperature, we calculate the fraction of the dust mass 
in several ranges of temperatures: above 50K (hot), between 25 
and 50K, between 15 and 25K and a cold 10K component for the 
two galaxies NGC 1705 and Haro 11. We estimate the fraction 
of dust between two temperatures as the mass of dust exposed to 
a radiation field such that the large silicates, which are at thermal 
equilibrium, have temperatures between these two temperatures. 
The energy absorbed by a dust grain is given by : 



/ 

Jo 



labs 00 



(4) 



where Q a b s is the absorption efficiency and U the mean intensity 
of the radiation field. 

The energy emitted by a dust grain is given by : 



Y em oc I v^B v (T eq )dv oc T*q P 
Jo 



(5) 



where B y is the Planck function and T eq the equilibrium temper- 
ature. 

As T a b s =T em , we can equalize these two expressions. We thus 
obtain a relation of proportionnality between the radiation inten- 
sity U and the equilibrium temperature T eq . Grains are assumed 
to possess an emissivity J3=2\ 



4+£ / T 



(6) 



Note that in this equation, the equilibrium temperature is nor- 
malised to the equilibrium dust temperature of the Galaxy of 

17.5K ( |Boulanger et al.|[T996l ). 

Finally, from the prescription of |Dale et al.| ( |200l] ) (Eq. [TJ, 
we can derive the heated dust mass associated to each radiation 
intensity U and thus to each equilibrium temperature. The frac- 
tion of the total mass exposed to radiation between temperatures 
Tl and T2 is given by: 



fmass([Tl, T2]) 



r 6(l-or) 



r 6(l-a) 
1 



U, 



l-a 



u 



l-a 



r 6(l-or) 



r 6(l-a) 



(7) 
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25 
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21 
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20 
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45 
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52 
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12 
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36 
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41 
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12 
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38 
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46 
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1.69 xlO 5 


2 


6.6 xlO 6 


2.02 xlO 7 


14 
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Table 8. Dust masses derived from SED models using different shapes of ISRFs (§ 6.4): a synthesised ISRF of a non-ionising cluster 



(created after an instantane ous burst, Salpeter IMF) and the synthesised ISRFs of 4 low metallicity galaxies (Galli ano et al.||2QQ3 
|2Q05l|Madden et al.[[2Q06] ). 

* The models used for Haro 11 and NGC 1705 include an independent, very cold dust component at 10K with an emissivity J3=l. 
M warm dust is the dust mass without including the 10K dust mass (but with cold dust included in the fit) and Md ust is the total dust 
mass of the galaxy, 10K dust mass included. 





NGC 1705 * 


Haro 11 * 


Mrk 1089 


Tgg (max) (K) 


88.2 


125 


146 


T e q (min) (K) 


10 


10 


16.9 


f m ^([50K;150K]) 


0.1 % 


0.4% 


0.15 % 


f m ([25K;50 K] ) 


7.2% 


28.9 % 


9.5 % 


U ( [15 K ; 25 K] ) 


8.3 % 


0% 


90.4 % 


W(10K) 


84.4 % 


70.4 % 


0% 



Table 9. Minimum and maximum equilibrium temperature and distribution of the dust mass with dust temperature ranges for the 
galaxies NGC 1705, Haro 1 1 and Mrk 1089. 

* The models used for Haro 1 1 and NGC 1705 include an independent, very cold dust component at 10K with an emissivity J3=l. 



NGC 1705 * Haro 1 1 * Mrk 1089 



Lr/i? - L[2ium,noOfim] (L Q ) 
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5.8 x 10 7 


1.7 x 10 11 


3.9 x 10 10 




(2) 


4.9 x 10 7 


1.4 x 10 11 


3.8 x 10 10 




(3) 


5.9 x 10 7 


1.6 x 10 11 


4.3 x 10 10 


L[3^m,50^m] (Ad) 




1.9 x 10 7 


1.2 x 10 11 


1.4 x 10 10 


£[50jum,100jum] (Ad) 




2.2 x 10 7 


4.1 x 10 10 


1.3 x 10 10 


L[l00/um,U00ium] (Lq) 




1.6 x 10 7 


9.6 x 10 9 


1.1 x 10 10 


LfiR = £[8j"m,1100yum] (A)) 




5.1 x 10 7 


1.6 x 10 11 


3.7 x 10 10 


L Ha (erg s" 1 ) 




6.34 x 10 38 (c) 


3.2 x 10 42 M 


4.57 x 10 41 (a) 


Size 




0.5' x 0.5' 


0.25' x 0.25' 


0.31' xO.12' 


SFR (M yr" 1 kpc" 2 ) 


(4) 


1.6 x 10" 3 


5.1 x 10" 2 


5 x 10" 2 




(5) 


9.4 x 10" 4 




5.8 x 10" 2 




(6) 


5.2 x 10" 4 


2.4 x 10" 2 


1.4 x 10" 2 




(7) 


5.7 x 10" 4 


3 x 10" 2 


1.8 x 10" 2 



Table 10. Luminosities, size and SFR of NGC 1705, Haro 11 and Mrk 1089. 

- Total inf rared luminosities (L TIR ) are estimated: 1) by integrating our SED models with a submm constraint from 3 to 1 100 //m, 2) 
from the |Dale & Helou| ( |2002| ) formula, 3) from the praine & Li ( |2007| ) formula. We also give the distribution of the luminosities 
within different wavelength windows. 

- Ha luminosities (L Ha ) are given in a) Canno n et al. ( |2006| ), b) [Bergvall et al. ( 2000]) and c) |Iglesias-Paramo & Vilchez]dl997|). 

- Star formation rates (SFR) are estimated: 4) using the relation of Kennicutt (1998 ) based on L FIR , 5) using the relation of |Zhu et aL 
(2008) based on v L y [24 /mi]. (This last relation is only valid for galaxies with a L TIR < 10 11 L and can not be used for Haro 1 1), 
6) using the relation of Calzetti ( 2007 ) based on the 24 yum luminosity (L24) alone, 7) using the relation of Calzetti ( 2007 ) based on 
L 24 and L Ha , . 

* The models used for Haro 1 1 and NGC 1705 include an independent, very cold dust component at 10K with an emissivity J3=l. 
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where U m / n and \J max are the boundaries of the intensity range 
obtained from the SED modelling. 

For Haro 11 and NGC 1705, the fraction of dust at 10K is 
given by the ratio between the dust mass of the cold dust compo- 
nent and the total dust mass of the galaxy . The results are sum- 
marized in Table [9] The cool phase (< 25 K) constitues the major 
part of the dust in our galaxies, at least 70% for Haro 1 1 and 
up to 90% for Mrk 1089 and NGC 1705. The SED of Haro 11 
peaks at usually short wavelengths (36 /mi), with a significant 
fraction (~ 30%) of the dust mass at a temperature > 25K while 
this same dust mass fraction does not exceed 10% for NGC 1705 
and Mrk 1089. Thus their global SEDs reflect the different lev- 
els of SF activity (and/or morphologies) and the consequences 
on the dust heating. 

6.6. Total Infrared luminosity 

We calculate the total IR luminosity (L^) for our galaxies by 
integrating our modeled SEDs from 3 /mi to 1100 /mi and com- 
pare our L TIR values with prescriptions in the literature using IR 
broadbands. For example, Dale & Helou (2002 ) have made phe- 
nomenological SED models to derive L TIR based on IRAS ISO 
and some 850 /mi data and provide a recipe to estimate the 3 //m 
to 1 100 /mi luminosity using the Spitzer MIPS bands: 



L TIR = 1.6 L 24 + 0.77 L 70 + 1.35 L 160 



(8) 

In the same fashion, |Draine & Li (2007 ) expand this L TIR 
relation to include the IRAC 8 //m flux. 

L TIR = 0.95 L 8 + 1.15 L 24 + L 70 + L m (9) 

We estimate the L TIR of our galaxies using these two formu- 
las and compare to the L TIR we determine from our SED model. 
For the galaxy Mrk 1089, the 160 /mi observation is not available 
but the flux density at 160 /mi can be extrapolated from our SED 
models. We estimate this flux density to be 5334 (± 266) mJy 
when using the submm constraint. For Haro 11 and NGC 1705, 
the Lti R was calculated from the revised models which include 
a cold dust component with an emissivity j3= 1 . Results are pre- 
sented in Table 10 For the UM 311 system, we derive, from our 
SED model with a 160 /mi constraint equal to our upper limit, a 
L TIR which should not exceed 1.4 x 10 9 L . 



The values given by the |Dale & H elou ( 2002 ) formula for 
the L TIR seem to be lower for our galaxi es than the L TIR der ived 
directly from our SED models while the praine & Li| ( [2007] ) for- 
mula leads to closer result s. To know if o ur overestimation of 
L TIR compared to the Dale & Helou (2002) formula is a system- 
atic effect observed in other galaxies, we derive the L TIR of the 
galaxies of Galliano et al. (2008) by integrating their SEDs from 
3 /mi t o 1100 /mi and compa re t he results to the L TIR obtained 
by both|Dale & Helou|([2002| ) and |Draine <OI| ( |20U7l ) relations. 
Engelbra cht et al. ( 2008) also calculate the L TIR of their galaxies 
from their SED modelling. We calc ulated the L TIR der ived from 
the |Dale & Helou| ( [2002] ) and the |Draine & U\^00l\ relations 
for their galaxies (from the fluxes available in the paper) and add 
these galaxies to the sample. The ratio between these different 
estimates are plotted in Fig. [6] as a function of metallicity. We 
clearly obs erve a systema tic underestimation of the L TIR when 
using the Dale & Helou ( |2002| ) r elation, with a shift of about 
20% and a scatter of ~ 20%. The |Draine & U\^2mi\ formula 
seems to better estimate the L TIR of the galaxies to within ~ 
20%. For the lowest metallicity galaxies used in this comparison 
(12+log(0/H)<7.5), we note a potentially systematic underesti- 
mation but larger samples would be required to investigate this 
point. 



Table 10 also presents the distribution of the IR luminosity 
according to wavelength windows: from 3 to 50 /mi, from 50 to 
100 /mi and from 100 to 1 100 /mi. For Mrk 1089 and NGC 1705, 
the luminosity is roughly evenly distributed over the wavelength 
windows while for Haro 11, 70% of the L TIR comes out at wave- 
lengths shorter than 50 /mi. While the longer wavelengths ac- 
count for the major fraction of the dust mass, not more than 6% 
of the L TIR comes out from wavelength superior to 100 /mi. 

6. 7. Star Formation Rates 

To investigate how these galaxies behave with respect to the 
Schmidt law, which describes the tight relationship between 
global SFR and gas density, originally formulated for normal 
spirals, we derive the st ar formation rate s (SFRs) for our galax- 
ies using the relation of |Kennicutt| ([ 1 998j) : 



SFR (Af yr~ l ) = 4.5 x 10 44 L FIR (erg s~ l ) 



(10) 



where L FIR is the IR luminosity obtained integrating our SED 
models between 8 and 1000 /mi. 



We comp are these values with the SFRs obtained by |Zhu 
|et al.| ( |2008| ) for star forming galaxies, formula which is only 
valid when h TIR < 10 11 L© and thus can not be used for the LIRG 
Haro 11: 



n rip — i \ vL y [24//m] 

SFR 24 , m (M Q yr ) = — - 



(11) 



We finally use the relations of Calzetti ( 2007) to derive SFR 
estimates from the 24 /mi luminosity L24 and the Ha luminosity 

(L#ar)- 



SFR (Af yr~ l ) = 1.24 x 10" 38 [L 24 (erg ^ 1 )] a88 
SFR (Af yr~ l ) = 5.3 x 10" 42 [L Ha + 0.031 L 24 ] 



(12) 
(13) 



The SFR estimates are presented on Table 10 The values 
given by the Kennicutt (1998) relation are higher than the one 
obtaine d with the |Calzetti| ( |2007| ) formula but compare to the 
ones of |Zhu et al.| ( |2008| ). 

We investigate the behaviour of our galaxies with respect to 
the Schmidt law ( |Kennicutt||1998|) We plot the SFR estimated 
with the relation of Kennicutt (1998) as a function of the gas den- 
sity for our three galaxies. We overplot a large sample of spirals 
and starburst galaxies presented in ( |Kennicutt||1998| ); (Fig. [7]). 
The two galaxies Mrk 1089 and NGC 1705 seem to follow the 
Schmidt law. For Haro 1 1, we take the large upper limit of the HI 
+ H 2 mass (~ 10 9 M Q - |Bergvall etaL 2000 ) into account in the 
gas mass estimate. This galaxy does not fall close to the Schmidt 
law. Its location in Fig.[7]might be due to the fact that it is a dwarf 
galaxy, which can explain why the gas seems to be consumed at 
a much higher rate than in normal spiral galaxies. Undergoing a 
merger, Haro 1 1 should also have a high star-formation rate with 
a lower gas co nsumpti on time than less l uminous galaxies (see 
|Kennicutt et al.[[T994l|Tacconi et al.||2006| for studies on normal 
nearby and submillimiter galaxies respectively). Nevertheless, 
|Kennicutt|(l998| ) suggests that the SFR in LIRGs do follow the 
Schmidt Law dependance on gas surface density, which means 
that the relation linking the SFR and h FIR of starburst systems 
should not differ significantly from normal disk galaxies. In con- 
clusion, the location of Haro 11 as an outlier confirms that some 
gas is missing and that the current HI and CO observations do 
not lead to an amount of gas coherent with what is usually ex- 
pected. The fact that Haro 1 1 falls far from the relation favors 
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Fig. 6. a) Ratio bet ween the total i nfrare d luminosty (h TIR ) derived by integ rating directly the SED m odels from 3 to 1100 fim and 
Ltir given by the |Dale & H elou (2002) formula. The squares indicate the Engelbrach t et al.| ([2008) galaxies whereas the circles 
indicate the Gallia no et al. ( 2008| ) sample. Stars represent our three galaxies Haro 11, Mrk 1089 and NGC 1705 in metallicity 
increasing order . b) Same as a) but performed with L TIR given by the|Draine & Li|([2007]) formula. 



the cold dust explanation over the change in the emissivity of the 
grains for this galaxy. To follow the Schmidt law, Haro 1 1 would 
require 10 10 M© of gas, which would imply an order of magni- 
tude more than that suggested by CO and HI measurements and 



which is coherent with the "missing" gas mass required to obtain 
a D/G expected by the chemical evolution models (see § 6.3). 
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Fig. 7. Positions of Haro 11, NGC 1705 and Mrk 1089 in the Schmidt- Kennicutt diagram. The triangles represent normal spirals 
and squares represent starburst galaxies (see Kennicutt, 1998 , for details on the galaxies represented by triangles). Red stars indicate 
our sources. The position of Haro 1 1 in this diagram accounts for the HI+H 2 upper limit of total gas mass (~ 10 9 M o ) of the galaxy. 



7. Conclusion and summary 

The quantification of the dust mass of a galaxy aids our under- 
standing of its evolution and star formation history. Larger dust 
masses are sometimes found in low metallicity galaxies when 
using submm constraints in the SED modelling. In this context, 
submm observations are clearly necessary to lead a more com- 
plete description of the distribution and properties of dust. 

We focused our paper on the dust modelling of four low 
metallicity galaxies observed with APEX/ LABOCA. 

In this paper: 

1 . We present the first images of four dwarf galaxies carried out 
with the APEX/LABOCA instrument observing at 870 fim. 

2. We construct the SEDs with Spitzer IRAC and MIPS bands 
as well as the IRS spectra for Haro 1 1 . We apply our SED 
model and determine the dust properties of these galaxies. 

3. We find that the mass of PAHs accounts for 0.08 to 0.8 % of 
the total dust mass of the galaxies, which is a factor of 5 to 
50 lower than that of the Galaxy. 

4. To investigate the influence of the submm constraints on the 
interpretation, we test the effect on the SED model results 
when submm 870 yum observations are taken into account 
and compare with SED models not taking into account the 
870 yum flux, but with observational constraints at wave- 
lengths only as long as 160 yum. We find that the use of 
submm observational constraints always leads to an increase 
of the total dust mass derived for our low metallicity galax- 
ies. 

5. We choose to include an additional component to account 
for the excess submm emission of NGC 1705 and Haro 11. 
A cold dust component (~ 10K) with a J3 emissivity index 



of 1 substantially improves the fit. We find at least 70% of 
the total dust mass residing in a cold (~ 10K) dust compo- 
nent for these two galaxies. We note that describing a cold 
component of (3=2 does not give very different ;f 2 values, but 
would give unrealistically larger D/G. Our results however 
do not rule out the hypothesis of a change in dust emissivity 
as a function of wavelength proposed in recent studies (e.g. 



Pupae et aLl|2003||Meny et al.[|2007) . 

6. While Haro 11 has a substantial (~ 70%) cold dust com- 
ponent, it also harbours a significant fraction of dust mass 
(30%) in a warmer dust component (> 25K). The SED peaks 
at unusually short wavelengts (36 /mi), highlighting the im- 
portance of the warm dust. 

7. We determine the D/G for Mrk 1089, NGC 1705 and Haro 1 1 
to be 1.9 x 10" 3 , 4.1 x 10" 3 and 0.2, respectively. For 
Mrk 1089 and NGC 1705, these D/G are consistent with cur- 
rent chemical evolution models. On the contrary, Haro 11 
has an excessively high D/G considering the upper limits de- 
tected in HI and CO. Haro 1 1 also falls far from the Schmidt 
law, perhaps due to the observed deficit of gas in this galaxy. 
This could suggest the presence of a large amount of molec- 
ular gas. 10 times more molecular gas, compared to that de- 
duced from CO measurements, may be present but not nec- 
essarily traced by CO observations. 

8. From our SED models, we determine the total infrared lu- 
minosity of our galaxy sample to range from 5.8 x 10 7 for 
NGC 1705 to 1.7 x 10 11 for the LIRG Haro 11. These values 
of Ltir are sy stematically higher than those obtained using 
the Dale & Helou| (|2002| ) formula but compare better to the 
|Draine & Li| { |2007} formula. While -90% of the dust mass 
is residing in the FIR to submm regime, not more than 6% 
of the total IR luminosity in Haro 1 1 emerges from the FIR 
to submm (100 to 1100 yum), while most of the luminosity 
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(70%) emerges in the NIR to MIR (3 jim to 50 jim) win- 
dow. This is in contrast to Mrk 1089 and NGC 1705 which 
distribute their luminosities more equally in these two wave- 
length windows. 

Better observational coverage of the Rayleigh Jeans side 
of the SEDs should help us to disentangle possible scenarios 
to explain the excess detected in some galaxies in the submm: 
dominant diffuse ISM dust, modifications of the dust optical 
properties at submm wavelengths, very cold dust component etc. 
The Herschel guaranteed time key program of 50 dwarf galaxies 
(PI: S. Madden), which covers a wide range of metallicity values, 
will provide a broader coverage of wavelengths between 60 jim 
to 600 yum and will enable us to better sample the warm and 
cold dust. The observations will especially enable us to study 
the slope of the Rayleigh Jeans side of the SED and to learn 
if an independant very cold dust component does indeed exist. 
A following paper will discuss the systematic increase in the 
dust mass estimate with or without submm constraints and the 
implications for the D/G ratios for a broader sample of galaxies 
and study how it could possibly be influenced by metallicity. 

Acknowledgements. We thank Xander Tielens, Ant Jones and Eli Dwek for 
very stimulating discussion which improved the content of this paper. We 
would also like to thank Gerhardt R. Meurer for his high quality HI map 
of NGC 1705. This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between the Max- 
Planck-Institut fur Radioastronomie, the European Southern Observatory, and 
the Onsala Space Observatory. 



References 

Aloisi, A. et al. 1998, in European Southern Observatory 
Astrophysics Symposia, Vol. 55, European Southern 
Observatory Astrophysics Symposia, ed. W. Freudling & 
R. N. Hook, 154-+ 

Arendt, R. G. et al. 1998, ApJ, 508, 74 

Bendo, G. et al. 2006, ApJ, 652, 283 

Bendo, G. J. et al. 2008, MNRAS, 389, 629 

Bergvall, N., Masegosa, J., Ostlin, G., & Cernicharo, J. 2000, 
A&A, 359, 41 

Bergvall, N. & Ostlin, G. 2002, A&A, 390, 891 

Bernard, J. R et al. 1996, in American Institute of Physics 
Conference Series, Vol. 348, American Institute of Physics 
Conference Series, ed. E. Dwek, 105-111 

Bernard, J.-P. et al. 2008, AJ, 136, 919 

Bottner, C, Klein, U., & Heithausen, A. 2003, A&A, 408, 493 
Boulanger, F. et al. 1996, A&A, 312, 256 
Calzetti, D. 2007, Nuovo Cimento B Serie, 122, 971 
Cannon, J. M. et al. 2006, ApJ, 647, 293 
Chen, H.-W., Prochaska, J. X., & Bloom, J. S. 2007, ApJ, 668, 
384 

Chen, H.-W. et al. 2009, ApJ, 691, 152 
Condon, J. J. et al. 1998, AJ, 115, 1693 
Conti, P. S. 1991, ApJ, 377, 115 

Conti, P. S., Leitherer, C, & Vacca, W. D. 1996, ApJ, 461, L87+ 

Dale, D. A. & Helou, G. 2002, ApJ, 576, 159 

Dale, D. A., Helou, G., Contursi, A., Silbermann, N. A., & 

Kolhatkar, S. 2001, ApJ, 549, 215 
Dale, D. A. et al. 2005, ApJ, 633, 857 
— . 2007, ApJ, 655, 863 
Draine, B. T. & Li, A. 2007, ApJ, 657, 810 
Draine, B. T. et al. 2007, ApJ, 663, 866 

Dultzin-Hacyan, D., Masegosa, J., & Moles, M. 1990, A&A, 
238,28 

Dumke, M., Krause, M., & Wielebinski, R. 2004, A&A, 414, 
475 

Dupac, X. et al. 2003, A&A, 404, Lll 
Engelbracht, C. W. et al. 2005, ApJ, 628, L29 
— . 2007, RASP, 119, 994 
— . 2008, ApJ, 678, 804 

Ferland, G. J. 1996, Hazy, A Brief Introduction to Cloudy 90, 

ed. G. J. Ferland 
Fioc, M. & Rocca-Volmerange, B. 1997, A&A, 326, 950 
Galliano, F., Dwek, E., & Chanial, P. 2008, ApJ, 672, 214 
Galliano, F., Madden, S. C, Jones, A. P., Wilson, C. D., & 

Bernard, J.-P. 2005, A&A, 434, 867 
Galliano, F., Madden, S. C, Jones, A. P., Wilson, C. D., Bernard, 

J.-P., & Le Peintre, F. 2003, A&A, 407, 159 
Gil de Paz, A., Madore, B. F., & Pevunova, O. 2003, ApJS, 147, 

29 

Gordon, D. & Gottesman, S. T. 1981, AJ, 86, 161 

Gordon, K. D. et al. 2005, PASP, 117, 503 

— . 2007, PASP, 119, 1019 

— . 2008, ApJ, 682, 336 

Gordon, M. A. 1995, A&A, 301, 853 

Greve, A., Becker, R., Johansson, L. E. B., & McKeith, C. D. 

1996, A&A, 312, 391 
Grimes, J. P., Heckman, T., Strickland, D., Dixon, W. V, 

Sembach, K., Overzier, R., Hoopes, C, Aloisi, A., & Ptak, 

A. 2007, ApJ, 668, 891 
Guhathakurta, P. & Draine, B. T. 1989, ApJ, 345, 230 
Guseva, N. G., Izotov, Y. I., & Thuan, T. X. 1998, Kinematics 

and Physics of Celestial Bodies, 14, 41 



22 



M. Galametz et al: The SED of dwarf galaxies using LABOCA 



Heckman, T. M. & Leitherer, C. 1997, AJ, 114, 69 

Heisler, C. A. & Vader, J. P. 1995, AJ, 110, 87 

Helou, G. et al. 2004, in Bulletin of the American Astronomical 

Society, Vol. 36, Bulletin of the American Astronomical 

Society, 701-+ 
Hickson, P. 1982, ApJ, 255, 382 

Hirashita, H., Kaneda, H., Onaka, T., & Suzuki, T. 2008, PASJ, 
60, 477 

Hoefner, S. 2009, Cosmic Dust - Near and far, Th. Henning, 
E. Griin, J. Steinacker (eds.) ASP Conf Ser., in press 
(arXiv:0903.5280) 

Hopkins, A. M., Schulte-Ladbeck, R. E., & Drozdovsky, I. O. 

2002, AJ, 124, 862 

Houck, J. R. et al. 2004, ApJS, 154, 211 
Iglesias-Paramo, J. & Vilchez, J. M. 1997, ApJ, 479, 190 
James, A., Dunne, L., Eales, S., & Edmunds, M. G. 2002, 

MNRAS, 335, 753 
Johnson, K. E. & Conti, P. S. 2000, AJ, 119, 2146 
Johnson, K. E. et al. 2007, AJ, 134, 1522 
Kennicutt, Jr., R. C. 1998, ARA&A, 36, 189 
Kennicutt, Jr., R. C., Tamblyn, P., & Congdon, C. E. 1994, ApJ, 

435, 22 

Kennicutt, Jr., R. C. et al. 2003, PASP, 115, 928 

Kunth, D. & Schild, H. 1986, A&A, 169, 71 

Laor, A. & Draine, B. T. 1993, ApJ, 402, 441 

Lara-Lopez, M. A., Cepa, J., Bongiovanni, A., Castaneda, H., 

Perez Garcia, A. M., Fernandez Lorenzo, M., Povic, M., & 

Sanchez-Portal, M. 2009, A&A, 493, L5 
Leroy, A., Bolatto, A. D., Simon, J. D., & Blitz, L. 2005, ApJ, 

625, 763 

Li, A. & Draine, B. T. 2001, ApJ, 554, 778 

Lisenfeld, U. & Ferrara, A. 1998, ApJ, 496, 145 

Lisenfeld, U., Sievers, A., Israel, F., & Stil, J. 2001, Astrophysics 

and Space Science Supplement, 277, 105 
Madden, S. C. 2000, New Astronomy Review, 44, 249 
Madden, S. C. 2005, in American Institute of Physics 

Conference Series, Vol. 761, The Spectral Energy 

Distributions of Gas-Rich Galaxies: Confronting Models with 

Data, ed. C. C. Popescu & R. J. Tuffs, 223-+ 
Madden, S. C, Galliano, F., Jones, A. P., & Sauvage, M. 2006, 

A&A, 446, 877 

Madden, S. C, Poglitsch, A., Geis, N., Stacey, G. J., & Townes, 

C. H. 1997, ApJ, 483, 200 
Marleau, F. R. et al. 2006, ApJ, 646, 929 
Mathis, J. S., Mezger, P. G., & Panagia, N. 1983, A&A, 128, 212 
Mauch, T., Murphy, T., Buttery, H. J., Curran, J., Hunstead, 

R. W., Piestrzynski, B., Robertson, J. G., & Sadler, E. M. 

2003, MNRAS, 342, 1117 

Meier, D. S., Turner, J. L., Crosthwaite, L. P., & Beck, S. C. 

2001, AJ, 121,740 
Meny, C, Gromov, V, Boudet, N., Bernard, J.-R, Paradis, D., & 

Nayral, C. 2007, A&A, 468, 171 
Meurer, G. R., Freeman, K. C, Dopita, M. A., & Cacciari, C. 

1992, AJ, 103, 60 
Meurer, G. R., Staveley- Smith, L., & Killeen, N. E. B. 1998, 

MNRAS, 300, 705 
Meurer, G. R. et al. 2006, ApJS, 165, 307 
Moles, M., Marquez, I., Masegosa, J., del Olmo, A., Perea, J., & 

Arp, H. 1994, ApJ, 432, 135 
Moshir, M. & et al. 1990, in IRAS Faint Source Catalogue, ver- 
sion 2.0 (1990), 0-+ 
Munoz-Mateos et al. 2009, ApJ, 703, 1569 
O'Halloran, B., Madden, S. C, & Abel, N. P. 2008, ApJ, 681, 

1205 



O'Halloran, B., Satyapal, S., & Dudik, R. P. 2006, ApJ, 641, 795 
Ostlin, G., Amram, P., Masegosa, J., Bergvall, N., & Boulesteix, 

J. 1999, A&AS, 137,419 
Poglitsch, A. et al. 1995, ApJ, 454, 293 

Popescu, C. C. & Tuffs, R. J., eds. 2005, American Institute 
of Physics Conference Series, Vol. 761, The Spectral Energy 
Distributions of Gas-Rich Galaxies: Confronting Models with 
Data 

Reach, W. T. et al. 1995, ApJ, 451, 188 
— . 2005, PASP, 117, 978 

Relano, M., Lisenfeld, U., Perez-Gonzalez, P. G., Vilchez, J. M., 

& Battaner, E. 2007, ApJ, 667, L141 
Richer, M. G. et al. 2003, A&A, 397, 99 
Rouleau, F. & Martin, P. G. 1991, ApJ, 377, 526 
Rubin, D. et al. 2009, A&A, 494, 647 
Rubin, V. C. & Ford, Jr., W. K. 1983, ApJ, 271, 556 
Rubin, V. C, Hunter, D. A., & Ford, W. K. J. 1990, ApJ, 365, 86 
Sanders, D. B., Mazzarella, J. M., Kim, D.-C, Surace, J. A., & 

Soifer, B. T. 2003, AJ, 126, 1607 
Shetty, R., Kauffmann, J., Schnee, S., Goodman, A. A., & 

Ercolano, B. 2009, ApJ, 696, 2234 
Smoker, J. V, Davies, R. D., Axon, D. J., & Hummel, E. 2000, 

A&A, 361, 19 
Stansberry, J. A. et al. 2007, PASP, 119, 1038 
Tacconi, L. J. et al. 2006, ApJ, 640, 228 

Terlevich, R., Melnick, J., Masegosa, J., Moles, M., & Copetti, 

M. V. F. 1991, A&AS, 91,285 
Tosi, M. 2003, Ap&SS, 284, 651 
Walter, F. et al. 2007, ApJ, 661, 102 
Weingartner, J. C. & Draine, B. T. 2001, ApJ, 548, 296 
Williams, B. A., McMahon, P. M., & van Gorkom, J. H. 1991, 

AJ, 101, 1957 

Wu, H., Zhu, Y.-N., Cao, C, & Qin, B. 2007, ApJ, 668, 87 
Yun, M. S., Verdes-Montenegro, L., del Olmo, A., & Perea, J. 

1997, ApJ, 475, L21+ 
Zhu, Y.-N., Wu, H., Cao, C, & Li, H.-N. 2008, ApJ, 686, 155 
Zubko, V, Dwek, E., & Arendt, R. G. 2004, ApJS, 152, 211 



